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ABSTRACT 


In this thesis a profitability model is developed for 
an exploratory well drilling program. The model considers the 
number of wells drilled, the probability of success, the reserve 
size discovered by a successful well, the discounted profit per 
barrel for oil produced and the cost of an unsuccessful well or 
"dry hole''. Each of these parameters are analyzed and a prob- 
ability distribution developed for all but the dry hole costs. 
With this model and the respective probability distributions, 

a FORTRAN IV computer program is developed using the Monte 
Carlo technique of stochastic simulation. The program provides 
the best estimate of discounted profit for the venture, the 
risk of the venture as measured by the standard deviation (and 
the coefficient of variation) and the probability distribution 
of the discounted profit - all after any number of prespecified 
model executions. 

After developing a hypothetical set of base parameter 
data and executing the program, a reference set of values for 
the profitability and risk are obtained. Each parameter set 
is then varied holding the remaining parameters at their base 
values. Using this technique the effects of parameter changes 
on observed model results are assessed. 

Several different sensitivity studies are performed in 
this manner using the data thus developed. These studies 
include the effects of parameter changes on the expected profit 
and risk of the venture using a Significance Ratio technique 


V 








developed in the thesis. In addition, the convergence of the 
best estimate of discounted profit with the number of model 
executions, as well as the relationship between the number of 
wells drilled and the venture risk are analyzed. 

Also studied is the possible effects of periodicities 
in the combined distribution data with an increased number of 
model executions and the effects of changing a parameter prob- 
ability distribution on the best estimate of discounted profit 


for the venture. 
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CeAri eke. 
INTRODUCTION 


fea General 

A primary consideration of management in making invest- 
ment decisions is the eventual profitability of the proposed 
venture. There are several ways of measuring profitability 
such as Rate of Return, and Present Worth or Discounted Expect- 
ed Profit. The main purpose in calculating this profitability 
criterion is to summarize into a single measure the quantifi- 
able factors affecting the economic desirability of the project 
under question [25]. 

Profitability calculations have generally been based 
on the best estimates of the individual variables involved. 
However, these best estimates may or may not occur. The re- 
sultant single value of profitability based on these best es- 
timates gives no information on the chances of doing much 
better or much worse then the expected value. 

In his anticle, Risk Analysis in Capita ny csimicine 
[23], David B. Hertz demonstrates that in certain cases the 
conclusions reached on investment decisions may have been much 
mirrerent™a& the analysis provided notwonly a single procit— 
ability estimate but also the chances of that estimate being 
exceeded or, on the other hand, not being met. Hertz concludes 
that the application of probabilities will often yield entirely 
different and better decisions. Thus, by developing the pro- 


bability distribution of the profitability, management can 








make quantitative assessments of the risks involved in approving 
a particular investment. The chances of achieving a minimum 
desired return or of losing money can be determined so that 
management knows the size of the risk it is undertaking. 

In order to determine the profitability of the venture 
and the probability distribution of the profitability, a model 
must be established containing all the factors which bear on the 
problem. In some cases the model is of such a nature that it 
can be solved analytically. That is, the model consists of 
known mathematical procedures or algorithms, which when manip- 
ulated, either directly or iteratively, provide a solution. 

Many types of managerial problems, however, are so complex that 
neither direct nor iterative analytical procedures exist for 
their solution, or, if such methods do exist, they represent 
sueh great difficulties in application that their utilization 

1s impractical. In such cases we turn to the use of a technique 
called simulation. Richards and Greenlaw [42] discuss simula- 
tion as follows: 

Simulation involves designing and utilizing a model 

which replicates some aspect of the firm's operations... 

In general, simulation models may be distinguished from 

other operations research models in that they are non- 

analytical in nature, calling for an experimental app- 


roach, in which workable, satisficing solutions rather 
than optimal ones are usually provided. 


In this thesis an 01il exploration investment decision model will 
be developed incorporating several different variables, all of 


which bear on the profitability of an oil exploration venture. 








Each variable will have a different probability distribution 
making the model so complex that analytical solution is not 
practical. We will therefore turn to simulation techniques, 
namely the Monte Carlo technique (described in Chapter II), to 
provide a measure of profitability and a probability distribu- 
mion Of that profitability. 

Before discussing the organization of this thesis it 
may be well to clarify one point which seems to cause a great 
deal of confusion, that is, the distinction between uncertainty 
and risk. Contemporary literature, by and large, still maintains 
F.H. Knight's distinction between these two terms [29]. Risk 
is characterized by situation in which the outcome is not known 
with certainty, but where the probabilities of the outcomes are 
either known or can be estimated. Uncertainty, on the other 
hand, refers to situations for which probabilities of the out- 
comes cannot even be predicted in probabilistic terms. 

Concerning this distinction, Alfred Rappaport [39], 
writes: 

While this distinction may be a useful pedagogical 

device, there is serious doubt about its applicability 

to decision makers in an organizational context. Here 
the distinction is blurred by the fact that decision 
makers generally have some feelings about the probabil- 
ities of future events. Admittedly, these feelings 
range from a degree of confidence bordering on certainty 
to ill-defined feelings approaching mysticism. The 
critical point remains that these feelings affect the 
judgments made by the decision maker. 

Based on the reasons expressed by Rappaport the terms "uncert- 


ainty'' and "risk" will be used interchangeably throughout this 


thesis. 








1.2 Organization of the Thesis 


This thesis treats the sensitivity studies in profit- 
ability analysis in two parts. The first part, beginning with 
Chapter II, describes the Monte Carlo method of stochastic 
simulation, a technique used throughout the thesis. In Chapter 
III a general profitability model for an exploratory oil 
drilling venture is developed. Chapter IV discribes some 
general statistical "yardsticks" used to measure central ten- 
dency and variability of probability distributions. Beginning 
with Chapter V and continuing through Chapter VIII each varia- 
ble, or parameter, in the profitability model is analyzed and 
a probability distribution best describing the parameter is 
established. In addition, a hypothetical set of values for 
each parameter is established which will be used as a Base for 
the subsequent sensitivity studies. With the profitability 
model and the respective probability distributions established, 
a FORTRAN IV computer program is developed in Chapter IX using 
the Monte Carlo method to provide the best estimate of dis- 
counted profit for the venture, the probability distribution of 
the profit and the risk of the venture. Finally, in Chapter X, 
using the hypothetical set of values for each parameter pre- 
viously developed, the model is executed and a reference, or 
Base, set of values are obtained for the venture profitability 
and risk. 

The second part of the thesis deals with the actual 


sensitivity studies using the model, parameter probability 








distributions and Base values developed in previous chapters. 
Chapter XI describes the need for sensitivity analysis 
in general terms and sets the framework for the studies which 
follow. In Chapter XII the results of the sensitivity studies 
are reported and analyzed. Conclusions and recommendations 
resulting from these studies and analyses are contained in 


Chapters XIII and XIV, respectively. 





CHAPTER II 
STOCHASTIC SIMULATION - THE MONTE CARLO TECHNIQUE 


aol History 


Through the ages man has used chance processes. 

Theories of chance processes were known in ancient times. The 
evidence is in the casting of the astragal bones from the ankles 
of sheep and goats in the manner of dice to subdivide land 
»olldleaase This was called casting lots. However, it was not 
until the 17th century that probability theory became formalized. 
Interestingly, this mathematical development arose out of in- 
Miinies into games Or Chance. 

Another important development occurred in the 1940's 
during war work on the atomic bomb. Scientists were frustrated 
by being unable to solve some complex nonprobabilistic mathe- 
matical problems directly by analysis. By combining chance 
processes and probability theory, they solved the problems in- 
directly. To do this, they created a new problem stated in terms 
of stochastic models in such a way that the new problem was equ- 
ivalent to the complex nonprobabilistic problem. But the sto- 
chastic version represented chance processes that could actually 
P-mcoLriIedeout.sethey actually performed. the chances processes 
represented by the probabilistic, or stochastic, formulas in 
the new version of their problem. Their methods were equivalent 
to flipping coins or throwing dice, but much more efficient. 

The results of these experiments were then substituted for the 


stochastic expressions of the new version of the problem to give 


6 





answers to the original nonprobabilistic problem. These scien- 


tists gave the name Monte Carlo to their new technique [5][22]. 


tee the Momce (Gar wer Leehnaquc 

Monte Carlo techniques were born when probability theory 
and actual chance processes were used to solve problems that had 
no stochastic aspects at all. Since that time Monte Carlo tech- 
niques have been used to solve a multitude of nonstochastic (or 
deterministic) type problems [33]. 

When dealing with stochastic (or probabilistic) simula- 
tion it is necessary to have a means for dealing with those 
probabilistic variables, the values which assume a frequency 
distribution. Richard F. Barton [4] defines stochastic simu- 
lation as: 

.«. One in which differing outputcs tuiala tose | eean 

be obtained without changing the inputs (ignoring random 

numbers as inputs). Specifically, this means that 

identical parameters, starting conditions, and input 
time path values produce varying outputs trial to trial 
and run to run. 
One commonly used method for working with stochastic variables 
in simulation models is again the Monte Carlo technique. The 
technique is one which both has application to an almost endless 
variety of stochastic simulation problems and at the same time, 
is relatively simple to comprehend and easy to utilize. S.W. 
Hess and H.A. Quigley [24] define the Monte Carlo technique as: 
a sampling procedure whereby complicated expressions 
involving one or more probability distributions may be 
evaluated. In essence, it consists of simulation an 
experiment to determine some probabilistic property of 
a population of objects or events by the use of random 


sampling applied to the components of the objects or 
EvVenes. 








A parameter value is obtained by drawing randomly from 
the parameter's probability distribution. In a similar 
manner, selections are made for each parameter value 
from its respective distribution. This set of parameter 
values is then substituted into the model and the first 
sample value of the independent variable is computed. 
Subsequent values of the independent variable are obtain- 
ed by drawing additional sets of parameter values. 

The individual values of the independent variable so 

obtained are in approximation to its true probability 

distribution. The approximate distribution approaches 
the true distribution as the sample size increases. 

The Monte Carlo technique actually involves a two sta- 
ge process. In the first stage, stochastic relationships be- 
tween the variables are determined in order to construct the 
Simulation model. These relationships may be determined from 
Surveys, company data, industry data, executive experience and 
the like. At any rate, once the relationships are determined, 
probabilities of occurrence of frequency distributions are 
assigned for each variable. The frequency distributions are 
then represented by associating groups of numbers to each value 
of the distribution, the size of the group being governed by 
the relative frequency of the value. As an example, assume 
60% of the items in the distribution have a value A, 30% have 
a value B and 10% have a value C. These three values can be 
represented by numbers between 00 and 99. Thus the items 
having a value A would be represented by the numbers from 00 
through 59, those having values of B and C would be represented 
by numbers from 60 through 89 and 90 through 99, respectively. 


Once the model is developed, frequency distributions 


assigned and numbers representing these frequency distributions 








determined, we move to stage two which involves running the 
Monte Carlo process to obtain a variable value. We do this by 
first selecting a random number (or numbers) from tables if the 
Simulation is to be performed manually or by building a random 
number generator into a computerized simulation program. 

(With some computer systems, such as the CDC 6600, it is poss- 
ible to make use of a random number generating function in- 
corporated in the computer systen itself.) Once the random 
number is selected (or generated) it is matched to the corres- 
ponding number assigned to represent each value of the frequency 
distribution of the variable in question. This in turn estab- 
lishes the value of that variable. This process is repeated 

for each variable in the model, the variables combined and a 
Single answer is obtained. Stage two is then repeated as many 
times as is necessary to obtain a statistically reliable sample. 
The various stages of the Monte Carlo process have been schem- 
atically diagrammed by M.D. Richards and P.S.Greenlaw [41] 


and are shown in Figure 2.1. 
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CHARTER sii 
THE PROFITABILITY MODEL 


bet General 
As noted in Chapter I, the profitability of an explor- 

atory well drilling venture and the probability distribution 
of that profitability will be analyzed using a model which 
incorporates the several different factors which bear on the 
problem, namely: the number of wells drilled in the progran, 
the probability of a success, the size of field discovered by 
a successful well, the discounted profit per barrel of oil 
discovered and the cost of a dry hole, or an unsuccessful well. 
The model essentially measures the profitability of the venture 
on the basis of Net Present Worth technique whereby the expect- 
ed cash flows are discounted at a suitable discount rate to 
reflect the time value of money. In the model, and throughout 
this thesis, the following notation will be utilized: 

m= lotal discounted promt oneany one dnl) line vendeavon. 

E(P) = The best estimate (mean) of discounted profit for 

the total exploratory well drilling venture. 


Number of wells drilled in the drilling venture. 


ro 
i 


p = Probability of a successful well. 

x = Number of successful wells, where x can assume values 
iON BOm ae Onn 

Q = Reserve size of oil field discovered by a successful 


well (in million bbls). 


ikl 








IZ 


YQ = Reserve size of oil fields discovered in the venture 


(in million bbls). 


e = Discounted profit per barrel for oil produced 
(in s7bbine 
oy) = Cost for an junsuccesstul swells (on dry sole. 


Using this notation, the model for a single drilling endeavor 
might be expressed as follows: 
P=xQe- (n-x) oo (lr) 
In an exploratory well drilling venture it must be 
realized that the number of successful wells is a function of 
both the number of wells drilled, n, and the probability of 
suecess, p. Therefore: 
x = f(n, p) (2) 
Dercner, the size of reserve discovered in the Veneume, yao. 
is a function of the number of successful wells. 
2Q = f(x, Q (3) 
It will be assumed that discounted profit per barrel, e, 
and dry hole cost per unsuccessful well, both as absolute 
values, are independent in that they are not directly derived 


from the number of wells drilled, etc. Therefore: 


P = [=Q = fixe Q)| e -Jon -[ x = fe pit Ch (4) 


In Equation 4, which will be used as the profitability 


model in this thesis, it can be seen that the discounted profit 


for the well drilling venture is equal to the size of reserves 








Ls 


discovered in the venture (which is a function of the number 
of successful wells and the reserve discovered by each success- 
ful well) multiplied by the discounted profit per barrel for 
0il which will be produced from the reserves, less the number 
or unsuccessful wells (a function of the number of wells 
drilled and the probability of successes) multiplied by the 
cost of the unsuccessful well, or dry hole cost. 

In terms of best estimates, the best estimate of 


discounted profit for the venture, E(P), can be written, 
E(P) = E(x) [ £(Q) E(e) | -[n - EG] ¢, (5) 


The relationship expressed in Equation 5, has been used 
bpeemany tO estimate the profitability of a venture directly as 
follows: 

Based on all available information, "best guesses" are 
made of the probability of success, which, in this case, is 
the chance that at least one oil field will be found, the 
magnitude of the reserves which would be expected by a success- 
ful well, the profit per barrel for oil produced and the dry 
hole cost. It should be emphasized that normally only one 
value for each parameter is provided, that is the best estimate 
(or guess). These best estimates are then combined with the 
number of successes determined by 

x = np (6) 
That is, if we had an 18 well program and the probability of 
success was 1 in 9 then we could expect 2 successes. If we 


multiply the 2 successes by the best guess of reserve size per 
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well (say 5 million barrels) we would have a total reserve of 
10 million barrels. If we estimate the most likely discounted 
profit per barrel to be $.80 then we would expect 10 million 
bbls x $.80/bbl1 or $8 million. From this we would have to 
deduct dry hole costs which would be (18-2) or 16 wells at, 
say, $600,000 per well, or $9.6 million. The best estimate of 
discounted profit for the venture based on this technique is 
therefore - $1.6 million. In all probability the venture would 
not be started on the basis of this analysis. 

As discussed in Chapter I, what this analysis ignores 
is the range of values possible for each of the parameters. 
For example the best estimate for reserve size was 5 million 
barrels per success. That may be the best estimate, but what 
if the reserves actually exceed this best estimate by 25%? 

The best estimate of discounted profit would then have been 
$0.8 million and a potentially profitable venture may have 
been passed up. 

In this thesis the basic relationship expressed in 
Equation 4 will be utilized for the profitability analysis. 

In subsequent Chapters each of the parameters will be explored 
to determine their functional relationships and proper frequency 


aestributions. 








CHAPTER IV 


DESCRIPTIVE MEASURES OF PROBABTEI yy sDisIRicolieN: 


4.1 Measures of Central Tendency (Mean, Median and Mode) 


In analyzing each parameter in the profitability model 
as well as the resultant discounted profit, it is necessary to 
know not only its probability distribution but also some single 
measure typical of the parameter. This single measure will be 
used to describe the central tendency, or the general region of 
Gemtinuum in which the distribution lies. 

The most common measure of central tendency is called 
imenMean or Expected Value and is designated Fe (Uae x) Or Ec 
The expected value is a measure of location in the sense that 
it roughly locates a middle or "average'' value of the variable, 
x [14]. The mean is the sum of a set of values divided by the 
number of values summed. Expressed mathematically the mean of 


a set of n numbers, X1, X2, eo oh Gy LEE 


X41 aL X9 + XZ a XA ses + Xn (1) 


E(@QReu = = 


Oimewnere £; 1s the frequency of each observation of x (1.¢., 


Xys X2, Xz, etc.) designated as x; we can write: 


0 es xe 


E (x) =X = + (2) 


A second measure of central tendency is the Median. 
The median of a set of values is the middle value when the 


observations are arranged in either ascending or descending 


15 
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order. The median is sometimes defined as the value above and 
below which an equal number of values occur [31]. On a cumu- 
lative probability (or frequency) plot, the median corresponds 
to the value opposite 50% cumulative probability (or frequency). 
The third measume ot centrametendency. 15 tiem modcrmmmic 
mode is defined as the most frequently appearing value. It is 
elsoeneterred)» to as the "“mostelikely valle”. @1t shomlidebe noted 
that a distribution does not have to have a mode (i.e. the uni- 
form distribution). Conversely, a distribution may have more 
than one mode (hence the designated "bimodal" for a distribution 


With two modes). 


Why is there a difference between these three measures 
of central tendency? Lawrence Mann[32] describes this diff- 
erence as follows: 


This difference is due to the fact that the distribution 
1s somewhat skewed. Skewness is the degree to which a 
distribution departs from symmetry. If a distribution 
leans to the right then it is said to be skewed to the 
left (negative skewness); if it leans to the left, then 
it is said to be skewed to the right (positive skewness). 
See Figure 4-1]. 
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Skewed right (positive skew) Skewed left (negative skew) 


Figure 4-1 Illustration of Skewed Distribution 
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In this thesis we will analyze the parameters of the 
profitability model in terms of their mean, median and mode, 
as appropriate. In addition, the expected value and median 
for each case, or parameter set, will be measured. The expected 
value (mean) will be considered the best estimate of discounted 
profit for each parameter set and will be used, along with the 
median, to Compare the various cases in the sensitivity studies 


which follow. 


4.2 Measures of Variability (Variance and Standard Deviation) 

The expected value of a random variable is an "average" 
value. However, this measure provides no information about the 
variability of the values of the variable. A measure of this 
variability, or the "spread" or "dispersion" of the values of 
the random variable, adds another dimension to the description 
of the random variable. 

The most useful quantification of uncertainty associ- 
ated with a future event is the standard deviation, o. It 
measures the extent of the spread between individual possible 
outcomes of the event and the expected value of the event. 
Another measure of variability is provided by the variance, 
found by squaring the standard deviation. 


For a discrete random variable, x, the variance is: 


e a2 
; ig (x; -x) 


Z=E£ (x-x)? PE (3) 


2 n-l 


Xx 
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where x = the mean value of x 

X; = any value On 

n = number of values of x 

i; = thequencyaet oc: 
The denominator in Equation 2 is replaced by n-1 when the 
statistics are carried out on a sample size n taken from a 
large population - when trying to draw conclusions concerning 
the large population. Division by n-1 in Equation 3 makes 


the sample variance, s A an "unbiased" estimate of the pop- 


x 
. : 2 
ulation variance hy 


If x is a continuous variable, the variance is: 


oo 
: ae 
a2 = B (x-x)" = f (x,-%)” g(x) dx (4) 
- 00 
Evaluating the usefulness of variance as a measure of dispersion, 
Samuel Goldberg [15] remarks: 
One difficulty with the variance is that it does not 
measure dispersion in the same units as the values of 
x. Thus, if x has dollar values, then E(x) is a certain 
number of dollars, but since the variance is the mean 
squane deviation, Vare(aieis Measia cde scaolerancs 
squared. It is in order to have a measure of disper- 
Sion in the same units as the values of x that we 
define the standard deviation as the square root of 
the variance. 
Hence, variance can be adjusted for this disadvantage by taking 
its square root, thus defining the standard deviation, o. This 
adjustment makes the standard deviation directly comparable 
with values of the random variable. The standard deviation 


takes into account, as does the expected value, not only the 


extreme values but all intermediate values, logically relating 
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uncertainty to the expected value [45]. 

In this thesis the standard deviation for the Base and 
all other parameter sets is computed by first calculating the 
variance using a relationship derived from Equation 3 and then 


taking its square root as follows: 


(5) 





meso Coefficient of Variation 

A knowledge of the expected value and standard deviation 
gives much information about a random variable, even though one 
does not have any knowledge of the individual values of a 
random variable. However, standard deviation is a measure of 
absolute uncertainty, that is, it measures the actual amount 
of variation present in the random variable; and it is depen- 
dent on the scale of measurement. There is a need for additional 
important information to facilitate comparison of variation of 
several random variables [44]. 

The measure which relates uncertainty to the expected 
value of the random variable is the coefficient of variation, 


cv, the formula for which is as follows: 


0 
CNe = 


(6) 


si |ye 


where cv denotes the coefficient of variation. Use of this 
measure of relative uncertainty makes the variability figures 


directly comparable. It describes the degree of uncertainty 
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or relative risk of the individual expectations about the ex- 


pected value of the random variable [38]. 


Weston and Brigham [47] describe the use of the coeff- 


icient of variation in the following passage: 


... Use the coefficient of variation to compare the 
riskness of alternative investments whenever the stand- 
ard deviation would be misleading. In general, the 
standard deviation is appropriate when percentage rates 
of return are considered, while the coefficient of 
Variation is appropriate if dollar returns are consid- 
ered. The reason for this distinction is that rates of 
return are already "standardized" for size of invest- 
ment, but some form of standardization is necessary 
when dollar returns are concerned. 


In the sensitivity studies which follow, the risk 
involved with each Case, or parameter set, will be measured 
in terms of its absolute amount (standard deviation) and in 
terms of its relative value (coefficient of variation). These 
values will then be related to the standard deviation and 


coefficient of variation of the Base set of parameters indicat- 


ing the change in risk associated with each parameter change. 








CHAPTER V 


PROBAB Til i ySOr SSveccess 


5.1 General 

In assigning the probability of occurrances of future 
events one normally starts by examining past experiences with 
Similar events. 

The petroleum industry has been very methodical yin 
recording drilling statistics for years. These statistics will 
provide the decision maker with the necessary data to compute 
probability of success, or, as they are called in the industry 
the national success ratios. These national success ratios are 
helpful but may be misleading. They include, for example, all 
wells that produced any amount of oil or gas, regardless of 
whether the amount was small or large, wells drilled in different 
formations, shallow and deep wells, and wells drilled with 
various kinds of skills. 

While the national success ratios provide one set of 
past experiences they may be of little help in assuming prob- 
tioeerties for a particular venture in a specific location. 

For that particular venture local experiences rather than 
national experience might be more helpful. But even with local 
eayerience there are still problems in makingwthe Statiseics 
selective so as to have relevance to the venture at hand. For 
instance, should all wells in all geological provinces be in- 
Cluded? Or only those that are producing from the target 


horizon being sought in the venture under consideration? 


‘asl 








Le 


These are difficult questions to answer and may be complicated 
by the fact that 1f we wanted to exclude all wells not producing 
from a specific geological province a breakdown of data per- 
mitting this exclusion might not be available. 

It should be noted that while local and national ratios 
of successes do provide guides, the assignment of probability 
of success for a particular venture is largely subjective, and 
in some cases, intuitive. But, no matter how the decision is 
arrived at, a probability of success for the venture must be 
stated if a profitability analysis is to be made. Before 
examining how the probability of success will be used in the 
profitability analysis, let us examine some concepts associated 
with probabilities as they relate to business decisions. 

For a long time only interpretation of probabilities 
has been the classical long-run, relative frequency arguement, 
1.€., repetition of the event over and over under identical 
conditions. For the majority of business problems, therefore, 
this concept has been almost useless as a practical guide for 
action. Most business decisions concern singular, or only 
occasional, events that may never be repeated - particularly 
under identical conditions. Therefore, choosing an act because 
it offers the best "long-run average" fits only a limited number 
of business decisions. Yet the "long-run" philosophy persists. 
During his research, Grayson [18] reported the following 


comments by operators: 
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If you drill enough holes, the law of averages will 
take care of you. If you keep playing the averages 
the percentages will hit sooner of later. 
When he questioned them as to what they meant by the words 
"sooner or later'', "average'', etc. most operators could not 
be specific. 

If we say that the probability of success is 1 in 7 
does this mean that if we drill 7 wells we will have one 
success? The answer is no! It means that if an event,such as 
drilling a well, is repeated over and over, the frequency of 
suecess will tend to approach the 1 in 7 success ratio. It is 
therefore important to remember that the number of successes 
is dependent on the number of wells drilled as well as the 
probability of success. If we consider the 1 in 7 ratio an 
average success ratio then, as Grayson writes [20] 

an operator who can drill only a small number of 
wells may expect large variations from the "average". 

The implications for a man with limited funds are clear. 

If he had unlimited funds, or even a very large amount 


of funds, he might be able to survive these variations - 
early losses or long "'runs'" of failures - while waiting 
for the "average to come in". But with limited funds, 
an operator may be ruined before the "average" is 
realized. It would be similar in example to the man 
who drowned while crossing a stream with an "average" 
depth of two feet. 
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Sezethe Binomial Ditetri bution 

An event that has a certain probability of occurrence, 
Dp, an a definite number of trials, n, is said to be binomial, 
mrscributed. 

When we drill an exploratory well for oil, one of two 
outcomes is possible. We can have a successful well, which as 
we will see in the next section is defined as a well which 
discovers a field of 1 million barrels or greater, or we can 
have an unsuccessful well, or a'dry hole'. If we let p equal 
the probability of success than l-p, or q , 1s the probability 
of a failure. When two exploratory wells are drilled three 
things can happen: both can be successful, one can be success- 
ful and one dry, or both can be dry holes. The chance that 
both are successful is pxp or p, that one is successful and 
Iieary 1S p X q + p X q or Zpq and that both are dry is q x q 
or q*. Using this same logic, if three exploratory wells are 
drilled we can have 3 successes (p>), 2 successes and one dry 


hole (3p2q), 1 success and 2 dry holes (3pq) or 3 dry holes 


(q>). Writing these terms in equation format we have: 


p> + 3p¢q + 3pq* + q> = (p+q)? (1) 


Thus, if n trials are made each of which results in either 
“swecess'' or failure”, if the probability of success on each 
trial is p and of failure on each trial is q=l-p, and if the 


trials are independent, then the probabilities of all possible 
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results of such a trial are given by the various terms of the 


binomial expansion. 


wel = -2 2 | : 
(ptq)P = p™ + np 7g + BU-Lyn-2, Tae eee 


Where the fourth term in Equation 2 is the general term and 
describes the probability of x successful exploratory wells 


in n independent trials. Hence, the binomial distribution: 


= n! ope nleg 
WOO = Srgmexyr Pod (3) 
The term a GE OH is often written as le (the number of com- 
binations). Equation 3 can therefore be expressed as: 
u(X) Seaeeep<q: (4) 


Again it is emphasized that Equation 3 provides the 
probability of x successes inn trials only if the probability 
Smmecuccess, p, 1S Constant from trial to trial and the trials 
are independent. 

To illustrate the use of the binomial distribution in 
predicting successful events let us examine the following 
example. 

Suppose we are planning an exploratory well drilling 
program of 20 wells and we estimate the probability of success, 
p, to be 1 in 7 or say 15%. Then we can also say that the 
beebability of failure, q is l-p, or 85%. Using Equations. 
we can write the probability of having x success in 20 trials 


aS.. 
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w(x) = Mose La BE) On (5) 


To determine the probability of 0 successes we would set 


x=0, and, 


w(O) = o9Co(.15)9(.85) °° 9 

_ 20! 0 20 
wO) = pADL (159°C. 85) 
alo) e fase 
w(0) = .0388 or 3.88% 


Likewise, the probability of 1 success (x=1) can be written: 


WO) = appar 159785) 


4) = 24329-x 1018 peti 

W = 2. oso. ae, ‘ P 
1.2165 x 10 

w(1) = 20(.15)(.0456) 

w(l1) = .1368 or 13.68 


A similar analysis could be performed for each success until a 
total of 20 were calculated. The terms could however, be read 
directly from tables [43]. For this example we would have the 


following probability for each possible outcome: 
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TABLE 5-1 
PROBABILITY OF SUCCESS - BINOMIAL DISTRIBUTION 


Number of 
Successes (x) w(x) in percent 


Ge 
1P32015 
22.99 
Lae 
VSeyZd 
10. 25 
54 
1.60 
. 46 
gl 
09 
. 00 


OOOO 


bo pes 
- FE OWWONKAWPWNFO 
[> 


20 0.00 


5.3 Simulation of the Binomial Distribution 

In the Monte Carlo studies to follow we will want to 
Simulate the binomial distribution in order to determine the 
number of successful wells for each drilling venture. For 
values of p near 50% and for large n it can be shown that the 
Binomial distribution approaches the Normal distribution [12]. 


The mean value, x and standard deviation, o, are given by: 


O°? 


vas 
Lf) 


o = MP (6) 


o = ynpq (7) 
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If xX, corresponding to a cumulative frequency, W(x) = 50%, 
whaeh 1S a Shorthand writing for W(x = x; ) = 50%, and ga e 
corresponding to a cumulative frequency, W(x) = 16% are plotted 


on normal probability paper and a straight line drawn between 
them, the probabilities associated with all other values of x 
could be read directly from the graph. Alternatively, the 
Normal distribution function 

1 ( (x-x,)4 


a : 
Oy /27 2c. (8) 


w(x) = 


could be solved analytically to obtain the value of w(x) for 
each success, xX. 

If the value of p varies substantially from 50% and if 
the value of n is relatively small can the normal distribution 
be used to simulate the binomial distribution? To answer this 
question we will continue the analysis of the 20 well drilling 
venture that had a probability of success of 15%. From Table 
9-1 a cumulative probability, W(x), can be calculated such that 


W(x £ xj) is the chance that the number of successes will be 


Xs or smaller. 





hws 
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TABEE $2 


CUMULATIVE PROBABILITY OF SUCCESS - EINOMEALT SDT ote eeu) 


Number of Successes (x5) w(x; ) in Percent (W(x =o aa Percent 
0 S.c8 Ga rells 
il 15205 Se 
2 22093 40.49 
3 24.28 ae thy 
4 sig Bl 92.98 
5 10.28 932216 
6 4.54 97.80 
7 oe 99.40 
8 0.46 99.86 
9 O11 99.97 

10 0.03 100.00 
11 0.00 LOGZU0 
20 0.00 100.00 


The values of (Wx £& x; ) corresponding to values of x 
on Table 5-2 were plotted on normal probability paper (Figure 
5-1) as circular points. If these values, in fact approached 
a normal distribution then these points would fall along a 
Straight line passing through xX, at (Wx =o) = 50% and a 


— 


at W(x =x,) = 84%, 


From Equation (6) 


X, = np 
xv = 2 Orel 5) 
x, = 3 
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From Equation (7) 


o = Y npq 

6 = Vv (AO ero 
O = y a2 soo 

o = 1.60 


Therefore, x, + o = 4.60 


O 


Peomaight line 1s drawn through these two points on Figure 5-1. 
It will be noted trom Figure 5-1 that the binomial 
@istribution points lie off the normal distribution line: It 


Can therefore be seen that for small values of p {in the range 


of 153) and for relatively small values of n, the use of the 


normal distribution would not be suitable to simulate the bi- 


iomeal distribution for the entire range of x;. 


With this in mind the binomial distribution for this 


study can be simulated indirectly by assigning, for each value 


of x, a series of consecutive numbers equal to the values of 


w(x) listed on Table 5-1. 


SIMULATION 


Number of Successes (x) 


— 


COW ONAN MN HP WN FO 


VRB IE oro 


OF BINOMIAL DISTRIBUTION 


0388 
1368 
VAs) 
2428 
Sal 
1028 
0454 
0160 
0046 
Oe ar 
0003 


w(x) x 1000 


0000 
0388 
150 
4049 
6477 
8298 
9326 
9780 
9940 
9986 
Shel el I) 


RANDOM NUMBERS 


Ose7 
ois 
4048 
6476 
Sage 
9223 
SS 
99359 
9985 
9996 
J 99 


Random Number Ranges 
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If we now generate a four digit random number from some random 
number generator this number can be equated to a specific value 
of x. See Table 5-3. As an example, if a random number of 
9395 were generated this would be equivalent to 6 successes. 
bikewise, if a number 0035 were generated this would equate to 
QO successes, or all dry holes. With enough trials the binomial 


distribution can be simulated in this manner. 


5.4 Coefficient of Variation 


From Equations 6 and 7 thescoefficient of variation can 


be derived as follows: 





eye 22 = = Siede (9) 
oe np 
Since q = (lp), Equation 9 can be wraeten: 
cvs — = 
From which we obtain: 
1 1- 
Clea 
‘na Pp (10) 


Using the data for number of wells drilled and prob- 
apbalaty of success previously noted, it is possible to calculate 


the Coefficient of Variation using Equation 10 as follows: 


SS as 


ih 
CV = Tae v5.60) 
20 lS 





st 


CV = 0.53 iE) 








Figure 5-1 
The Normal §& Binomial Distributions 


Thi re 
---ft-—} fit} -—— 4nane 





ese +4 ERG HHIHHGEEAL OBEEGEEEEEGELG oY 
BeBe a 4 gas) f- 
=| ESESS Basen suseeenses ab Srtaanp H+ Be 
Sela Shab ams to Saas a 
~ SoS Se o2aeea seee= Seta SS SS SSS SS SS geeoe Bio === =a = ts Se oe a BS Sa SS2 2am 


- ST TTT TTT eee 


Hse NEL ANEL EEE RE REELECTED 
SUOURGGHURERREREREEORSOCOREEEEREOOES” 400R08 REREROEES 
oe He aH oe HH 4H aH aoeee SeESE? se” "Seenuae BBa a 


OO 


Bees aia Se” ..aseeen See eae 

Bap pt tt BDOED SRSA PVURSSSESRADH’ saeaaDes oe lteta BOBROUB 

a Bea See Eeeseeceee ae el ee ey DRAPE SXSSS SSS SSE SS SRY 4245 SOOT Sa BS BeMRBeaaoeeG 
ee | ang Sage enT Eee FEEEee Saale tere S080 ST SRE0 FC NCORr sa SCOR RREOOS MA RAR see 
| Seas eunes seeeceeuacseeecceseesessnseeeeeeecers BRCSP 4B BERTA CARES 
Pies aes ap pel SRR USRSRRREEOoo SRR eo SECC eee Fy ves oo exe {2aee aeeoe 
0 Ct SES OSS NN OS O05 SSS EE PERS SSeS ESS SASS Cee 


[2] 
a an ALL LUTOUENETL VOT CTT 


99 


oy 





BORER TAOTRILAOMRRS IS Hoth oH BEBE f FECES 
Sep epee ones PERE EERE GH BEROE PRBRARE Dh eRe eeo 
EI Dpal ee oe eeeees Bee eS eeene reas ‘aha tH SB te EREREReS EAORR SER 
Sake Seah eeecee Gia Ht R287 ABReR ES RGE RRAREOREURERE IER 
i as | | a a | Sete ae tH SEH oe ¥ ADURSRUSS PEE EE \_} —-—----- + OGBbOSSREe8 

fr ee 0 GM se tlt EE 25.) (2980 50n0R 

a Ow Fates 

Bn a a ed) hl hl 

DRED) GARGS eeSSaR GDeap Seshens 
RRR S25 BSS SORES SSB aaa 








10 


90 





= 


HH 





aoe ‘ 
+} os a: iat aa © 
SUBUS SEEBSERSSU EST 22s GHEE SEEEREEESG GEER ZaG88 ee ss 
ROS SERBEeSeReCrS2Hh coe PERERA ESEE 
BESR’ 4828 oo BEER EEE Ee 
ERR BeES |---| +--+ -—-- —---—_-++—-} | © 
sneeee BSe200 +h as CGE BER GUOODOOERSSGBbEaGR BBSOoO Se0E02:0:0BSGH8 |b 
he | td T0982 Case Cee eee Seas SSSSSREEOO Vocammee I. aS eea 
: ses SEP EHS eRe 





20 


30 





40 


specs EgEEEEEEEEEEEEEEEEEEHEEEEEEEFEEEEECEC 


Zz _ BRS eS | 4 BES Shep SERB SaSaa tt Bae @ wn 
oon ‘oiaieseoos Ooo oa0n oases Sau8 Benoa 
7 a SBRERaea cena 6 BOR GSaDS0GGE0 042 ae@aee GESSBSaS00 OSeeea | 
mY r. p— + | +4 an BBe GEE SaaS SaSeeaR faeae Beet pit +4 BSGO0 SeSR000G088 
SG Seelae BBOoLt aia nan Hee He nap —~-—- aS (ea a HGWE DEBS aRCan iS 
GERD Ht at a ae BSS Sooesheean et SSURa Ky 
28 SSNw’,.2ooe iatstaiaeate teletatapatay ame GeooC ideale Sq sseaeeas SOROS B8E8 BSR ROReS 
an SG’ 46600) SSRRODGoeo eeweo ORR ee a es eae ca | 0) ea 4i—-+- eee eae ele re 
Be’ 46BER SHRRRGEBRES DARED Peete Gemma BEREBERBEE SEALS =s OSE 
GGa5 saacs a? JES BESSERE RSS BRDOE BSRSeS OOO SSE RBS Gee SEBREo0 peices Sa BOnoe| eS 
F d2SR ORO ER Ses CLES PLoS CARSALES +H G@W@@cnenhoGa eeeaa Slate 
SES oo5 Peo tH BSSR8) RRS GRRE Hetet t+ 4 HH Hy ttt tt . HH aa8 








60 








BRBELS BERR RG SREB RRR SRS 

peel SSGSRSGG0 GUUG0 COUGe BOUSe Ree 
°o J BERS GERGEBREEGME BS & o 
[2] a z e—- eae SE Bom N 








Se auae ana sart SETETSzs#otsovezstzerssosetverasezesasssas 

Bae eee HR Soe SSS eo See GEGB = SCR SSF Esso sees = 

a 88 PaO ee ea oe eae HEES HABIBI ATSEO GSU aRea8 
G0 FSESa VI OeE RECA LARA BERRA COREE ee 

SSR TPERRRER RAE SEE ERE EeER re 


ann 
Sep daeet eee teee onsa000050 iitiineettonnan HERE RESRRRE RM ai TH a TORanouanata 
oy a EEE HTH Onn TTT TT 


95 





98 





Pree FEEEEEEEEEEEEEEEE RECEEEEEEE EEE EEE EEEEEEE 
asEane Sis esssseons fosattasti fost aatoniiiee 


99 
1 











GBB 
OO 
Ng 
. Sela BetUS RRames Me ees HGR RASS BARR Bw 4 tH aan as 
Be PURE RSSEES BRoeeooose HH He SBSaneab alatatalal CB 
See See SGas tueeercussceeerceeee feaseseus sat tHe 
rSsgauue sEBBASu@UECH OO SOmaaE CHOEORGRUGRUWUHUUE ABeg 
“AEOSEBA ED ata AP RRRHE OLAROES Bai elaeeL ela AC 0 GEE 
“ Trent EUHEEUT aT RTEG HATH ATEL TTT RET ATTA ALATA RTE : 
sane AEG EEEEEaeL TT LE 


0.05 0,1 


BEreuerad teiaiuah PETE 
: PE lik AN, 2 
oa) 1 2 3 4 5 6 7 8 9” 
xX 


i) 


W(x < X 








CHAP TERS 
OTD KRECERY Oli oleae 


It goes without saying that the economic success or 
failure of a drilling venture may well be determined by the 
Size of the field discovered as a result of a successful well. 
Of course the size of the field is unknown prior to drilling. 
It is therefore essential, if we are going to attempt to es- 
tabilsh the relative profitability of the venture, to determine 
two items. First, whether there is a general mathematical ex- 
pression which can be used to characterize the probability 
distribution of field sizes in the geological province. In 
other words, do the field sizes vary in such a way that their 
probability of occurrence can be predicted? Secondly, if a 
general field size distribution function can be established as 
a model, what are the boundary conditions we should establish 
for our particular venture? Each of these questions will now 


be explored. 


6.1 Field Size Distribution 

At first glance it would seem quite unlikely that a 
general mathematical model of field size distributions could be 
established - considering the complex process of oil reservoir 
formation, or, if by chance a model could be established for 
one geological province it would have general application in 
another province. 

some ae the earliest work in this area was reported by 
J.J. Arps and T.G. Roberts [3] in 1958 after a study of reserve 

SS 
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data from the Denver-Julesburg Basin. In their study, a 5.7 
million acre sample was chosen on the east flank of the Basin 
(located in Nebraska, Wyoming and Colorado) and the reserve 
estimates of all 338 fields found to that date in the Lower 
Cretaceous sands were investigated. Arps and Roberts arranged 
the fields in nine groups according to their areal extent and 
estimated the ultimate primary recovery from each group. The 
limits of these groups were chosen in such a way that each 
group represented fields about twice as large as the preceeding 


one. The nine groups thus arranged are represented on Table 6-1. 


TABLE 6-1 


OIL RECOVERY DATA ON DENVER-JULESBURG BASIN FIELDS 


Range of Productive Number of Estimated Average 


Group Acres Fields Ultimate Recovery 

(M bbls) 
1 28 56 iLO) 
Z 23-57) aa 5 50 
5 oe = els 4] 16540 
4 PS) 12726 70 29 au, 
5 226 - 453 40 957.6 
6 A55.- 905 55 1 OV iGrez 
ii 905 - 1810 19 Bites Ln (6 
8 1810 - 3620 1 11,460.0 
9 3620 1 50.7 500 


When a frequency - density distribution of the data 
presented in Table 6-1 was plotted versus Average Ultimate 
Recovery on logarithmic-probability paper a straight line can 
be best fit through the resulting points. (See Figure 6-1). 
This demonstrates that the field sizes in the Denver-Julesburg 


Basin are Lognormally distributed. 
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While a Lognormal distribution best fits the data in 
the Denver-Julesburg Basin there is still no reason to general- 
1ze and say that all reserves follow a Lognormal frequency dis- 
tribution. If, however, the Lognormal distribution remains 
invariant under a wide variation in the definition of reserves 
and geological regions then a generalization can be made. 

To test this sinvariancy Gordon M. Kautman, (2/7) 4) inva 
extensive study reported in 1962, accumulated data from several 
different sources including: The 011 and Gas Conservation 
Board of Alberta Report to the Lieutenant-Governor-in-Council 
as summarized in The Financial Post; Oil and Gas Journal stat- 
istics on two states for 1946 and for 1960; an individual 
company's estimate of reserves in four basins; and the Arps - 
Roberts article on the Denver-Julesburg Basin. 

The data accumulated by Kaufman, covering fields in 
four coastal basins, fields in North and South Louisiana (both 
in 1946 and in 1959), fields in Oklahoma, etc. were plotted on 
sixteen (16) separate graphs on logarithmic-probability paper. 
Kaufman, in analyzing the results of these plots, writes [28] 


The graphs displaying this data tentatively suggest 
thar 


1. Within a given area the functional form of empirical 
histograms of reported field sizes does not change 
(a) over time, 

(b) with changes in the definition of 'reserves", 

(c) with changes in the minimum size of field 
reported, or 

(d) from one geographic area to another. 


2. The* functional form is the same for a well-defined 
geological basin such as the Denver-Julesburg as for 
an arbitary geographic area such as a state. 








oy 


Although each of /the sixteen (16) graphsm/ isvamplen 
of data from fields of varying geological ages and of 
varying length of production lives, the invariance in 
the functional forms is striking enough to warrant the 
conjecture that added data blocked on these two factors 
would not alter the conclusions stated above. 

A visual observation of these graphs reveals that most points 

generally fall along a straight line indicating a Lognormal 

frequency distribution of field size. 

In another analysis of oil field size distributions, 
Folkert Brons reports the result of a study of the predicted 
ultimate recovery of 75 Southern Louisiana Miocene Fields [10]. 
Pemonown On Figure 6-2, when the data is plotted on logarithmic. 
meobability paper the points fall along a straight line indi- 
cating a Lognormal distribution function for these fields. 

On the basis of the evidence presented in the studies 
cited above, of reservoirs with varying reserve sizes, geologic 
origin and physical characteristics, it will be assumed for 
further analysis that a Lognormal distribution function can 


Serve as a general mathematical model for reservoir size 


estributions. 


6.2 The Lognormal Distribution 


The lognormal distribution can be defined as the dis- 
tribution of a random veriable whose distribution obeys a 
normal distribution function. Thus, if x is lognormally dis- 


tributed then y = 1n x is normally distributed and, 


1 -(y-Yo)? 
wy) = —~— exp} ————— 
Oy V2T ; 20,7 7 (1) 
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If thewdistribution  tsedesitredmimeterms of xX rather 
than y a change in variables is performed which results in the 


probability density function for y expressed in terms of x as 





follows: 

- Z 2 

1 (Pix ve 

w(x) = CAD See (2) 
a, V2T “lone 
x 

where: x = independent variable 

o., = standard deviation 


yo 
Mee cane ae NO Dy 


The equation gives a lognormal distribution curve, 


with its characteristic positive skew, as shown below. 


ee 





Lognormal Distribution - Positive Skew 


Figure 6-3 
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Since y is normally distributed: 


y, (median = y , (mean) = y(mode), (3) 


but since x = e”%, x is not normally distributed and the values 


Ose x,, (median), x, (mean) and x(mode) must be calculated separ- 





ately [9]: 
The Median, Xn! 
y = In x 
y 
so Ge (4) 


x =e (4A) 


and, In x, = y, (4B) 


The Mean, x,: 


Yo + oy 
Se 2 (5) 
ay? 
= + 
ene Ale 2 We 5 (5A) 
The Mode, Xpoqge: 
_ one Z 
Seis oP (6) 


Oy eat a Oy, (6A) 
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purther, the Variance, o is given by the following expression: 


Xx 


Gane 
a4 = exp f{ Zt a4 (em - 15) 


oe 
insenerime Equation 5 an, Eiq@ualt 1 Orme sc Olea iin: 
2 
o* = x 4 (e%¥"-1) (8) 


X O 


From which the Coefficient of Variation, CV, may be expressed 


as: 
(9) 
0 Z 
Gy =e? ee 
O 
Ror Small values of Oy we can approximate Equation 9 as: 
CV = Oy (10) 


Om Boundry Condastaons 


If we assume that the field (or reservoir) sizes in 
any area are lognormally distributed and will fall along a 
Straight line when plotted on logarithmic-probability paper 
then the next problem is one of specifing the position of the 
line which would represent the reservoir size distribution in 
the specific location where the drilling venture is to be 
conducted. 

This line could be drawn if someone could predict the 
median value of the reserve size to be discovered by a (yet 
undrilled) successful well (i.e. the reserve size corresponding 
Hora SOS probability) and the standard deviation of the dis- 


terputrTon representing the logarithm of the field sizes. 
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Then, by inserting these values in Equations (4B), (5A) and 
(6A) previously noted for median, mean and mode a straight line 
could be drawn on log-probability paper. The problem is, how- 
ever, that the operator (or his agent) would probably not be 
able to specify either the median or the standard deviation. 
What he most likely would be able to specify is his estimate of 
the highest and lowest values of reserves he would expect with, 
say, 90% certainty. He would probably base these estimates on 
either his own experience or the experience of others who have 
drilled in the area. As an example, the operator may not be 
able to specify an exact reserve size, although he might say 
that he would expect the reserves to run from 5 million to 150 
million barrels with 90% confidence. If these two points were 
puoctedson log-probability paper - with the 5 million barrel 
figure plotted at the 5% probability point and the 150 million 
barrel figure plotted at the 95% probability point a straight 
line could be drawn between them which would represent the 
Lognormal distribution of all reserves in the area. 

With the line plotted, the median value of Q (or x in 
Equation 4B), corresponding to a 50% probability, and the value 
of Q median plus one standard deviation, corresponding to an 
84% probability, can be read directly. In addition, to use this 
model in our determination of the relative profitability of the 
venture using the Monte Carlo technique, it is essential that 
we be able to determine a reservoir size and corresponding value 
of probability (that the reservoir size will not be exceeded) 


for each point along the line. 








Oo 
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To illustrate this procedure, as well as to generate 
aaa, tor use in the Monte Carlo studies conductedwas part of 
this thesis, let us assume we are going to drill for oil on- 
shore in the Southern Louisiana area. By going to the approp- 
riate issue of the Oil and Gas Journal [48] we would find a 
listing of the sizes of all fields in the Southern Louisiana 
on-shore area. These fields could then be arranged in ascend- 
ing order of size (Q) and a relative frequency density, w, and 
cumulative frequency, W(Q < Q;) calculated. This tabular 
listing, along with the calculation of relative and cumulative 
frequency is shown as ACTUAL Values in Table 6-3. 

Having the actual values of cumulative frequency cal- 
culated, these values can be plotted on log-probability paper 
with the ordinate being the cumulative frequency and the ab- 
scissa being the field size. The triangular points on Figure 
6-4 represent a plot of the actual data for the Southern Louis- 
lana on-shore area. 

The operator can now use these points (or any other 
data he has available) as a guide in drawing a straight line 
on the log-probability paper beginning at some upper value of 
Q, corresponding to a 95% probability, and at some lower value 
@£)@, corresponding to a 5% probability. In this example the 
points chosen were 79 million barrels as the upper value and 
5 million barrels as the lower value. With this line drawn 
(Curve A, Figure 6-4) it can now be said that, for the purposes 
of this profitability analysis, 90% of all reserves discovered 


by a successful well will fall along the straight line, 5% of 





— = = 
> a 
ed 
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them will exceed the upper value and 5% will be less than the 
lower value. Theoretically then, there is a 5% possibility 
that a reservoir could be discovered with infinite reserves on 
one extreme and a 5% possibility of a reservoir with a 0 res- 
emvye on the other. Practically of course, this is not reason- 
able since 0 reserve means no oil and an unsuccessful well. 
If this straight line represents the reserve size distribution 
for a successful well then some lower limit must be established 
which would represent the reserve size below which the well in 
question would not be considered successful, even though some 
011 was discovered. In this example a lower limit of 1 million 
barrels was chosen. Therefore, we say to be considered succ- 
Gssitul a well must discover at least a 1 million barrel reserve. 
At the upper extremes it would be foolish to say that 
an infinite reservoir could be discovered even though this is 
possible by an extrapolation of Curve A, Figure 6-4, to a 100% 
Probability point. Realizing that this is not physically 
possible an arbitrary upper limit should be established. In 
this example it was observed that in the Southern Louisiana 
on-shore area the maximum field size is 137 million barrels 
(see Table 6-3, ACTUAL). On this basis an upper limit of 140 


Million barrels was chosen. 


6.4 The Lognormal Distribution of Field Sizes 


With the boundary conditions now established and Curve 
A,Figure 6-4 drawn, it is now possible to proceed with the 


calculation of the mean, median and mode as follows: 
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The Median, Xn: 


Q at 50% = 20.0 million barrels 
Len. Qnedian aan 


From Equation 4B 
es DS In aS 


ime 0 010 


<< 
O 
Il 


fi 2) 86 


X< 
T 


Yo (yam In x (at W< x = 84%) 


Mee - oy = In 46.0 
2.996 + Oy = 3.829 
Oy = 0.833 
a 
Oy = 0.694 
ae 
Oy Ono 
2 


The Mean, x,: 


From Equation 5A 


a2 
ae ar 
ln Xo = 2.996 + 0.347 
In > S2o05 
x = 28.3 million barrels 
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The Mode, peradic: 


From Equation 6A 


Ee 7 Z 
In orc Yo Oy 


iat Se sis 2.996 - 0.694 
ln aaa = 2.302 
de Me bs, LOO million barrels 


From Equation 9 


O 
oy = = = ody? 4 
oy V0. Sean 


CV ay 2 Syl 0 la 
CV = 1.00 


If we assign random numbers to represent W(Q < Q;) we 
can, from Curve A, Figure 6-4, manually determine the corres- 


ponding value of Q as follows: 
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TAB Ore 


FIELD SIZES - RANDOM NUMBERS 


WUC) Ss | Q 
Random No. in %) x 10°bb1) 

05 05 5.0 
06 06 ae 
07 07 ape 
93 93 69.0 
94 94 73.0 
95 95 80.0 


While the process of manually reading the graph, tab- 
ulating the values of Q corresponding to Random Variables 
which represent W(Q < Q;) and inputing these numbers as data 
in the Monte Carlo computer program is possible, the process 
1S quite laborious. To equip ourselves for an efficient sim- 
ulation system a Lognormally Distributed Variable Generator 
can be used to generate the values of Q on a high speed 
digital computer. 

Such a Generator was initially developed by McMillan 
and Gonzalez [34]. A version of the Generator, modified for 
use in this problem, incorporating the concepts of upper and 
lower limits, is included in Appendix A. Using this Generator 
it is necessary only to input the median Value of Q and the 


value of Q at W(Q < Q;) = 84%, both easily obtainable from 
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the plot of a straight line connecting the 5% and 95% points, 
and the upper and lower limits, as previously discussed. 

With this data inputed, the Generator will provide field sizes 
(Q) which are lognormally distributed within the boundary 
conditions established for the program. 

In order to test the reliability of the deviates 
provided by the Lognormal Generator, 1000 values of Q were 
generated using the program included in Appendix A and the 
boundary conditions previously discussed. The results of this 
test, along with a calculation of relative and cumulative 
frequency, is included on Table 6-3 under the heading of 
SIMULATION. For comparison purposes the values provided by 
the Generator were placed beside the corresponding actual 
values. Additionally, Curve B, Figure 6-4 was plotted with 
points generated by the Lognormal Generator. 

By comparing Curves A and B and the actual data points 
on Figure 6-4 it will be seen that excellent correlation 
exists except in the lower values of Q - but even there the 
difference is not considered significant. On this basis,the 
Lognormally Distributed Variable Generator will be used to 
Simulate lognormally distributed values of field size in this 


thesis. 





eG [9 S00 100 


G0 


NZ ol 


0S Ov o€ 


09 


02 


86 


66 


666 866 


66 66 





Figure 6-4 


Cumulative Dis 


tribution of Field Sizes 


Southern Louisiana On-Shore 


TTT iy Dy ) 


i {41 












= 
= | 
eo 4 
3 

es 
=| 


ees ees 


Prmessoces 


CEE RTEET ERIE Ti i EEA 


tb 


aT HLET ESE 5S SS Se i822} 


re 


a 
: 
He iaieea 


+ 
+ 
SEE EEE ETT EET 
+ 
a 
sans 
PLP Ee EEE EE EE 


+ 
+ 
+ 


GUTTER TERRE PE UE 


ba 


HEE TE 


a 

+ 

ie 
+ 


i 


¢ 
t 
+ 


teat 
a 
ett 
+ Wes 2X 


3. 
at | 
a @ 


2 
Bett eT 


—_—=- w 


003 
eae 
men 
eral 
igeceal 
ee 
Ea 
Rass 
oe 
Sane 
++ 
me 
eaae9 
= 
eo 
eo 
oes) 
at 
a 
of 
pits 
pease 


2 
= 


a 
fa 
HEE 
EL HEE HHT St 


. 


Pat Fpd PET 
SEE HIE 
ges: 
Sess 
» —4 


PSS =F LTT HT Fa) PT NNENT LEE EEE 
E 


sRetiii 


tL, | | 


ee fee) 
eee ease 
eee ssees epee sue: 
4 1 


HELE : 


eae 
STRUUREEEENG 


sal 


BeSSe c0vas 

BREE ELC iota Tye] IS Wee 
MESSSSC SIs a> Pel. | 
SSSESnO Teas? ia 
Pa ee | 


SSSR ERR...) com 


annee 


t 
' 


HECHT 

< HH 
mT 
an. Fr 
i) a) 

freq iH BaER 


fT ESREL ETE EE 
} 
| 
HE Es 


Nf 


a 


? 
4 


ESET Bae Bi PB SIzE Mh. 


> 
Q 
ct 
Cc 
seb) 
— 
PEER HT HEE TTT 


Se 
ca 
HE 


H 
He 


i 


a i 


biti! 








$bt6 666 


> OM 


(Fd 








20 


TABLE. 625 


SIZE OF FIEED - SOUTHERN LOUISTANA ON -SHnere 
(1970 Estimates) 


ACTUAL | SIMULATION 


Field Pred Re a iaee de) (C Winemaie ae Freq. Rel.sFreq: Gumeiacae 
Saez -Q 





fe 

; (£)  (w=px100)  (W <Q) | CE) (wemxl00) (WE Q) 
(x10°bb1) (3) (3) (3) (3) 
Hl 3 Bees OS 
2 lal eee Lae 
3 18 Late G2 
4 25 DiS Bs 
5 24 a il Bred 
6 2 a ostth 224 55 ree 11.4 
‘i 3 5 5 5.9 gl Se | 14 ll 
8 2 2 tl 8.35 34 eG Lees) 
9 0 ee 19 1.9 1oe4 
10 5 Sa 14.2 45 EPs) 2529 
11 4 anes 18.9 34 Sig tt Dye 5 
eZ 5 5.9 Zar. © las es. 29.8 
13 0 Dl ies 25 Zea Seal 
14 bs) 5.9 7 i Se 5526 
1S 7 SeeZ, 38.9 28 2.8 Giles Ie 
kG 3 Sao 42.4 5/ Sar uh 5 as 
1 4 ule AW |. 52 See WS 
18 s) 5.9 oo 0 Z5 es) 48.0 
19 Z eet So. 4 Z 0 Lao B05 
Zi 1 We 56.6 Cy Hed So aV4 
elk 1 2 Sih gee. all eel: SoS 
Ze 1 Ike Z 59.0 18 Hate Sel 
ZS 4 wy Oro... / ZZ peat 597.5 
24 iH als 64.9 20 20 (6) eS 
Z5 0 64.9 Al ieee ea, 
26 1 Wad Gier. 4 19 9 64.5 
fil 1 1 ee ‘(a AS lee ey Spats) 
28 5 52.5 jes: 12 i lots G0) 
29 il eee To) 13 sO 68.3 
30 il eee, 2 20 220 70S 
Si 1 i eae pal al 15 oes Tel ate 
Si oe 8 0.8 77-0 
33 pala ILS das: 73.9 
34 fe 18 te: Goal 
55 Z 2.4 TOS 6 O26 TOs 
36 1 ead fio 518, ILS) jes: TAs 
oy Poe 0 4 eee Fe .0 
So ore O 8 Q.8 78.8 
39 eo 7 0 F/ oS 
40 2 2G oO) al 6 6 oOo 
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TABLE 6-3 Continued 


ACTUAL SIMULATION 














Field Freq. RelFrequa) Cunmeced: Freq. Rel.Freq.  Gun@erec- 
Sze -Q £ i 
6 (f) (w=7x100) (W < Q) (CE) CWS a1 0105 (W < Q) 
feol0° bbl) % ) (% (%) (% ) 
4] 1 1h Ys 81.6 5 0.5 80.6 
42 Zz a i 84.0 8 0.8 81.4 
43 84.0 e) 0.9 BiZie-5 
44 84.0 5 Obes Bi7e..0 
45 1 eZ Sra 6 0.6 3, Oe 
46 852.2 8 0.8 84.0 
47 il jie 86.4 4 0.4 84.4 
48 86.4 9 0.9 oes 
49 86.4 4 0.4 85.7 
5) 86.4 2 iw 85.9 
SJ i 1.2 Seno y 0.9 86.8 
Se Be 5 0.5 Si 
53 1 fe 88.8 5 OES 87.8 
54 88.8 Z Oy 88.0 
5 88.8 6 Oe6 88.6 
56 88.8 5 Ome5 89.1 
oy 88.8 2 az SOeo 
58 88.8 i) 0.7 90.0 
59 1 ee 90.0 5 Des 90.3 
60 Z Dae 92.4 6 0.6 90.9 
61 92.4 6 0:6 91.5 
62 92.4 2 Oe 91.7 
63 92.4 7 Oiney. 92.4 
64 il 2 93.6 2 OreZ 92.6 
65 93.6 1 0.1 92.7 
66 93.6 5 O55 93.2 
67 93.6 1 Oe 93.3 
68 93.6 2 eZ 955 
69 93.6 4 0.4 93.9 
70 93:.°6 5 O35 94.2 
Gl. 93.6 3 0:..3 94.5 
ae 93.6 5 OES 94.8 
UG 93.6 94.8 
74 93.6 5 O35 95.1 
75 1 Ja 94.8 1 ore: 95.2 
76 94.8 2 O22 95.4 
a 94.8 1 O21 95.5 
78 94.8 il 0.1 95.6 
79 94.8 il Or Oi] 
80 94.8 95.7 
81 94.8 1 O.a8 95.8 
82 94.8 5 O25 96.1 
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TABBE 6-3 Continued 


ACTUAL SIMULATION 






Field Freq. Rel.Freq. Cum.Freq. | Freq. RKel.Pweq. 9 Cumm@iied. 
paze -Q £ £ 
: (f) (weixl00) (WSQ) | (£)  (we_x100) (WE Q) 
oo Pb) (3) (3) (4) (4) 
83 94.8 96.1 
8 4 94.8 1 Oat 96.2 
35 94.8 4 0.4 96.6 
86 94.8 96.6 
87 94.8 96.6 
88 94.8 2 Cee 96.8 
89 94.8 96.8 
90 94.8 1 Oa 96.9 
91 94.8 96.9 
92 94.8 96.9 
93 iL ee Tone i Ost OF 0 
94 96.0 3 Gass OFS 
95 1 12 Oia 4 0.4 Oe 
96 O72 2 Oe 97.9 
97 97.2 97.9 
98 ae. 97.9 
99 O72 97.9 
100 97.2 97.9 
io) il Oy 2 eZ Oar 
102 97.2 1 eat 98.2 
OS OiaZ 1 Ore 98.3 
104 Cie 98.3 
Os ey 1 Ol 98.4 
106 ied iL end 98.5 
nO 97.2 98.5 
108 i Meayz 98.4 Oia 
109 98.4 Ocee5 
110 1 | ee 99.6 98.5 
att 99.6 98.5 
ez 99.6 98.5 
113 99.6 il Orr 98.6 
114 99.6 il On 98.7 
iS 99.6 98.7 
116 99.6 98.7 
ae 99.6 1 Oa! 98.8 
118 99.6 98.8 
ALS) 99.6 il Oca 98.9 
120 99.6 98.9 
ZA. 99.6 98.9 
222 99.6 98.9 
AZ 3 99.6 5 Ons 99.2 
124 99.6 99.2 
es 99.6 i Oeil 99.3 
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TABDES6-3 SGontinued 


ACTUAL SIMULATION 


Field Freq. Rel.Freq. Cum-.Preq. | Freq.) Kell {Prec ee ace 
See = 10) 





; (f) (wegxl00) (WE Q) | (£) (we—xl00) (WE Q) 
ot) (8) (4) (4) (4) 
£26 215) 018 fle) & 
eZ / she) eee 
128 Ne) [6 Oo) fro 
eo 99 20 oh) 6 & 
130 oo OO 
131 ETRE S Ee S 
1b) 4 el se. o) eae) 
155 S18) ole eel & 
134 eee Ook 99.4 
135 SENS, elas 
1S hehe 99.4 
37 1 ie 100.0 eel 
SiS 100.0 1 Ore eae) 5 
es 9 POOe V rae) 
140 BOO 0 5 Deyo 100.0 
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CHAP TER yin 


DISG@UNTED PROFIT) PERSBARREE 


meee ceneral 

in a very real Sense, the wltimate proti Capnmercy oem 
exploratory well drilling venture, where oil is discovered, is 
dependent not only on the size of the new reserves but also 
(and perhaps more importantly) on how much profit the operator 
feeetereceive for each barrel of oil produced. The case can be 
Visualized where the operator in a successful exploratory drill- 
ing venture discovers vast potential oil reserves but the 
Sapemce Of production far exceeds the posted price he can 
jmeeeive fOr the o1l produced. In this case no matter how large 
meicmpotential reserve, the profitability of the venture is less 
mmeane0, Or 1f produced, the operator would suffer a loss. the 
cases of off-shore drilling in deep water and production on the 
North Slope of Alaska (the latter receiving wide publicity) 
illustrate this notion. It is therefore essential in a profit- 
meee analysis to estimate the expected profit per barrel on 
reserves discovered by a successful well. In practice, though, 
it is rare that the operator would commit himself to "the" 
expected profit per barrel for the venture. He could perhaps 
provide his estimates of the lowest, highest and most likely 
Values of profit per barrel basing these estimates on his 
experience in the area or on factors similar to those which 
permitted him to estimate the probability of success (see 


Chapter V). 
54 
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7.2 Discounted Profit Per Barrel Estimate 

In estimating the highest lowest and most likely values 
Se profit per barrel the operator would first have to provide 
his estimate of the highest, lowest and most likely posted 
Mmmmecemmor the oil. From this he would have to deduct his ex- 
Bemses including, royalties, severance taxes, operating costs, 
aeieriang and completion costs, etc. This would then provide 
his net profit per barrel for each or the three estimates. 
But, Since reserves cannot be produced all at once, the profit 
per barrel can not be figured on the basis of current sales 
prices. Production will be spread into the future (say 10 years) 
and future money is not as valuable as present day money. To 
meeemmt tor this financial fact-of-life the operator will have 
to convert his net profit per barrel calculated at today's 
prices into a discounted cash flow. In this discounting (or 
deferment) process a discounting factor is developed by calcu- 
lating the Present Worth of a $1 annuity (i.e. constant income 
or expenditure) received each year for a specified length of 


time, n, at a specified discounted rate, i, [16]. This is given 


wen foe AL) 
Loe 


where: P is the Present Worth 


by the expression: 


R is the amount of the annuity (taken as 


$1 in this development) 
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(lt an -] 1s the uniform series Present Worth 
ace 
-} 
1(1 1) or Present Value Factor. 
The Present Worth Factor can also be written as 


aa or Values of the Present Worth can be 


"7 a@+a 
i} 
obtained from Tables. If these values are then divided by the 
number of discounting periods, n, the deferment or discounting 
factor is obtained. 
To illustrate, the discount factor based on a discount 


rate of 10% for a period of 10 years is calculated as follows: 


Prom Equation. ; 


eles 1 yt 


le Sl 
i = 10% 
n = 10 years 


From Tables [17] the Present Worth Factor is 6.1444 


P = 1(6.1444) 
P = 6.1444 
Bascount factor = : Cs 
Discount factor = ote 
Discount factor = .6144 
Concerning this technique, Grayson [19] writes: 
Many operators, ... Object to thespresent iy aivcweonec pr: 


First, there are those who scoff at the idea of making 
any present value computations on exploratory ventures. 
There are already so many uncertainties, is their 
belief, that to refine the calculation any further is 
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like tilting at windmills. The answer to this has 

been stated several times: a present value calculation 

is an estimate, but an educated one that explicitly 

recognizes time as a realistic decision consideration. 
to ignore the effect of timemometo try tComeopom ie) mee 
in an implict way is to incur the danger that it may 

be overlooked entirely or its effect miscalculated. 

On the basis of the concepts discussed in this section, 
the highest, lowest and most likely values of profit per barrel 
can be computed. Table 7-1 illustrates this process and is 
based on the hypothetical expenses and prices for a drilling 
venture in the Southern Louisiana on-shore area. 

The values developed in Table 7-1 will be used as the 


base figures for the sensitivity studies contained in subsequent 


Chapters. 


me ewines triangular and Beta Distributions 


It was seen in Section 7.2 that while the operator 
would probably be unable to provide the exact profit per barrel 
of oil produced he could provide his estimates of the highest 
profit per barrel he would expect, the lowest profit per barrel 
and the profit per barrel he feels most likely to occur. It 
is assumed of course that these figures will be discounted to 
reflect the time value of money. Given this data, the question 
then becomes, what distribution function best suits this set of 
parameters? If all that is known is the upper limit, the lower 
limit and most likely value (mode) it is possible to describe 
this knowledge by means of either the triangular distribution 
minction or the beta distribution function. Of the two, the 


triangular function is the simplest and will be discussed first. 
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TABLE meek 


DISCOUNTED PROFIT PER BARREL - SOUTHERN LOUISIANA ON-SHORE 











MOST 

ITEM LIKELY HIGH LOW 

POSTED PRICE $3.25 3.50 3.00 
DEDUCTIONS: 

ROYALTY: 

@1/6 54 .50 

@1/8 44 
SEVERANCE TAXES (see Note 1) 24 AS ee: 
OPERATING COSTS (including lift 

cost, labor, lease maint., 

Bec...) .10 lb .09 

HOME OFFICE COSTS OR OVERHEAD 

(@ 5%) .16 .18 mals 
GEOPHYSICAL COSTS Oi 0 Fei 
DRILLING § COMPLETION .40 .40 .40 

(see Note 2) 

LEASE ACQUISITION (see Note 3) 02 .02 .02 
TOTAL DEDUCTIONS =¢ ee [Sle2  - 20 
NET PROFIT 1 ae 1.88 81 
DISCOUNTED PROFIT (based on $.88 $1.15 $.50 





discounting at 10% for 10 
years: discount factor = 
.6144) 


NOTE 1: [49] a. State and local severance taxes paid in 1968 
= $2550) 7 O00. 


D. Total Petroleum Praqurd preauction in Wooy— 
2,526,008 bbl/day (average). 


c. Average field price of oil (1969) = $3.33 


$233,071,000 7 
7,526,008 bbi/day x 365 = §-252/bb1 


Be i 





meme 2: [50] 
Pote 5S: a. 
[oye 
ome 
d. 


a. 


be 
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TABLEM?- 1 Contamued 


however, at $3.30/bbl 
Severance taxes = :422 = i255 4 
CeeSi) 
say 7.4% of posted price 


Assume production at 150 bbl/day. 


Estimated Cost of drilling and equipping wells 
in 1970 = $196,500. Allowing for inflation 

@ deeper wells assume cost of drilling = 
$217,500/well. 


Assume recovery of 217,500 in 10 years then 
recovery/year = $21,750. 


WAL Isle 


150 bbl/day x 365 = $-40/bb1 


Assume Lease Costs at $60/acre where a completion 


1S made. 


Assume well spacing at 40 acres; 


then Lease Costs $2400/well where well is completed. 


Assume lease on unproductive land to be equivalent 


to $2400/well; 
then total lease costs = $4,800/well. 


Assume recovery in 5 years with production rates 
150/bbl/well. 


$4800 
5 





= $960/well 


$960 


150 bbl/day x 365 ~ ‘91? 


say $.02/bbl 





_ 
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7.4 The Triangular Distribution Function 


The triangular distribution is defined by three values 
of the variable; i.e. the lowest valur of x = a, the highest 
value of x = c and the most likely value of x = b. The prob- 
ability distribution function w(x) can take three forms as 


shown below: 


w(b) 


w(x) 





aX 


Figure 7-1A Triangular Distribution - Positive Skew 


w(x) 


a b C 


Figure 7-1B Triangular Distribution- Symmetrical 
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w (b) 


w (X) 


i ee ee 


Figure 7-1C Triangular Distribution - Negative Skew 


C 
Since Hi w(x) dx = 1 (= the area of the triangle) for a 


Bomtinuous probability distribution: 


fw) dx = 1/2w(b)(c-a) = 1 (38 
from which 
2 


IS) aan (4) 


In unpubilshed lecture notes Folkert Brons [11] 
proves the following properties for the triangular 


distribution: 


ofa eet Te 
mean, X = = ae (5) 





mode, w(b) = z : ry (6) 


median, x > a (7) 


for positive skew (see Figure 7-1A) 
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x, = a + j{c-a) (b-a) (8) 


for negative skew (see Figure 7-1lc) 
Brons has also derived the cumulative distribution function 
[ll] and has demonstrated that: 
FOneeex = 5 
(x;-a)? 


W(x < x3) = (pra) Te-ay we 
or, xX. =a +\V W(x < eel lesa) (10) 


1 


For x < b 


ye 
a ay 
ON Xx; = ¢ - J W(x > Xx; ) (c-b) (c-a) (eZ) 
RuUrther , 
NOSES, | ile er CLS 


The variance for the triangular distribution is given by: 


2 a2 + b* + c4 - ab - ac- be 


= = [So a a (14) 


From which the Coefficient of Variation may be derived as 


moOlLLows: 


(ep 
< 
il 

i 
4! [bd 
il 
5 
-} 

—_ 

CO 
ia 
(@) 
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2 
a/ se [OB os ae - ab - ac - bec 
2 
CV = a t+ + C 
' es ae 
_ Va" + b’ + Cc” - ab - ac - be 
CY = 7 Sd Gar a) (15) 


foe oimulation of the Triangular Distribution 





Using the data for the highest, lowest and most likely 
values of discounted profit per barrel previously developed it 
is possible to simulate the cumulative triangular distribution 
of profit per barrel W(x) as follows: 


let a. = $250 


b = $.88 
ea SS 
x = specific values of expected discounted 


profit per barrel 


then for, $.50 < x; < $.88 
_ = .50 + W(x < x;)(.88-.50) (1.15 -.50) 
, = .50 + VW(x < x3) (.2470) (16) 


Bor, $l.l5 2 ee $88 


G - V W(x > x; ) (c-b) (c-a) 
as 1.15 - VW(x > eo nO oe eee: ae ete 
x. = 1.15 - VW(x > x.) (1755) Gwe 


ra 
it 


ra 
tl 


ra 
it 


ra 
tI 
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ehiial 


Wi xX >) 3 


W(x < $.88) = P-pES EG 


W(x £ $.88) 


il 
on) 
Co 
on) 


(18) 


Further, from Equation 15 the Coefficient of Variation may be 


ealculated as follows: 





cy = a? + bé = ae 2 Bley = Be = [oe 
7 1.414 (a + b +c) (4553) 
CV = 
564 
CV = 37580 
eyo = .157 


A series of random numbers can now be assigned to represent 
values of W(x < x3) seomeyil se OOO Aa a to .999999...... From 
Equation 16, if this random number was .585 or less the random 
number would be substituted into Equation 10, to obtain a 
Value of Xa. If the random number was .585 or greater 

W(x < x; ) = 1.0 - W(x > x; ) (19) 
would be used to obtain W(x > a) and then that value would 
be substituted into Equation 12, to obtain the corresponding 


value of x;. This process is illustrated in Table 7-2. 








DELUSTRAT ION: 

Random 

No. 
5 BHO. 00 (1. 
e011 .01 ( 
2 07 (2 
S58 2516 ts 
~59 (.59) 
. 60 (. O10) 
.99 (.99) 
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TABLE 2/22 


DISTRIBUTION 


PROFIT PER BARREL ASSUMING TRIANGULAR 


Maa we Uae Xx, = ) See 


i 
Witxes< x; ) W(x > x.) Vwex < Ke ciara) Vwex > Xa C7) 


00) 
99) 
98) 


42) 


41 
.40 


sel 


In this way the best estimate of 


2 Sit 
255 
Peis 
Hei) 

Bis), 

5 fhe 

| Peat Cal 


discounted profit per barrel, 


corresponding to a random number generated to represent a value 


of W(x < x;), may be simulated. 


For use in the Monte Carlo 


Simulation, a Triangular Distribution Variable Generator was 


developed and is included in Appendix B. 


In order to test the reliability of the deviates pro- 


vided by the Triangular Distribution Generator, 1000 values of 


the profit per barrel, e, were calculated using the program 


included in Appendix B and the values of a, b, and c previously 
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discussed. The results of this test are listed on Table 7-3 
under the heading of Triangular Distribution. The values of 
profit per barrel thus simulated were plotted on Figure 7-2 as 
Saveles. Also included on Figure 7-2, as triangular points, 
are values of profit per barrel calculated manually using the 
technique outlined on Table 7-2. It is apparent from Figure 
7-2 that there is excellent correlation between the manually 
calculated values and those calculated using the Triangular 
Distribution Generator, thus confirming the reliability of the 
Generator as a simulation tool. Dispite this excellent correl- 
ation one item should be recognized. By the very nature of the 
random number generator which generates values of numbers used 
to simulate W(x < x; ) a value of 100% for W(x < ey) can never 
be achieved. Instead the highest value of W(x < xj) we can 
achieve using the generator is .9999..... Thus, there is no 
possibility that the highest value of profit per barrel (e) 

can ever be reached. This shortcoming is considered minor 
however, considering the excellent correlation for virtually 


the entire range of deviates. 


7.6 The Beta Distribution Function 


In this derivation of the Beta distribution function 
Brons [13] notes that a disadvantage of the Triangular distri- 
bution is that its derivate ew Ge) becomes undetermined at xX = 
a, b and c, i.e. the only points where some information about 


x is available. This means that at the outer limits a and b, 








Figure 7-2 
Cumulative Distribution Curve - Triangular Distribution 
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wMixesuddenly drops to zero. This is depicted in Pieurcs7 3. 


w(x) 


dw (x) 


Xx 
—_—_—_—_—___> 


fmpure 7-3 - Triangular Distribution Derivative 


me also notes that a more gentle dropping off of the probabil- 
ity of occurence of x near the end points often better ex- 
presses reality - in that occurances in a real situation would 
normally taper off rather than end abruptly. This concept can 


be expressed by the Beta Distribution shown in Figure 7-4, 








a) 


w(x) 


a b c 


_— oe) S X 


Figure 7-4 Beta Distribution nine: nen 


M@mesBeta distribution, like the Triangular distribution, is 
defined by the three parameters a, b and c. 
The Beta Distribution is given by 


als: Fc CUES @ i 
ce) ol aay C2107) 


w(x) = 


It can be seen that the Beta distribution is a ratio 
of gamma functions where a and 8 are defined by the values of 
a, b, and c. A derivation of the relationships between q and 
8 and a, b and c (and of the Beta distribution itself) is be- 
yond the scope of this thesis and can be found in several 
excellent discussions of this distribution [13][19]. 

In an unpublished report [30], W.G. Lesso developed a 
computer program which generates deviates having a Beta dis- 
tribution. The program uses as inputs, the mean and variance 


of x (as expressed in a PERT type analysis): 
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(21) 


— 6 (22) 


From which the mean and variance for the Beta distri- 


bution may be obtained as follows: 


= 2 Gomer | 


x, ~ cra (25) 
Z 
i eee 
bo 2 
(c-a) (24) 


The Lesso developed program, modified for use in this study, 
is included in Appendix C. 

For comparision purposes the Beta Distribution Variable 
Generator was used to calculate 1000 values of the profit per 
barrel, e, using the same parameters of a,b and c as were used 
to Calculate the profit per barrel using the Triangular Dis- 
tribution Variable Generator. The values of profit per barrel 
thus calculated are listed in Table 7-3 under the heading of 
Beta Distribution. Cumulative frequency curves for the Beta 
and Triangular distributions, using the data contained in Table 
meeeis plotted on Figure 7-5. It can be seen from Figure 7-5 
that very little difference exists between the distributions 
except for values at the lower end of the curves. It should be 
noted, however, that the mean and standard deviation for the 
Beta distribution was computed using the values as expressed in 
a PERT type analysis (see Equations 21 and 22) rather than in 


-the form used in the Triangular distribution (see Equations 5 





oo 
= Be 
a re 
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ama 14). 


Because of its relative simplicity a Triangular 
distribution will be assumed for the variable, profit per 
@errel. In the sensitivity studies which follow however, the 
Beta distribution will be substituted for the Triangular dis- 
tribution and its effect on the venture profitability will be 


analyzed. 





TABLE 7-3 


UZ 


COMPUTER SIMULATED DEVIATES FOR DISCOUNTED PROFIT PER BARREL 


motit/bbl 
($) 


50 
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54 
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56 
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58 
ee) 
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ASSUMING TRIANGULAR AND BETA DISTRIBUTIONS 


parameters 


lowest profit per barrel: 


most likely profit per barrel: 


highest profit per barrel: 


TRIANGULAR DISTRIBUTION BETA DISTRIBUTION 
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TABLE 7-3 Continued 


TRIANGULAR DISTRIBUTION BETA DISTRICT Te 


. Relative Cumulative Relative Cumulative 
Profit/bbl | Frequency Frequency Frequency Frequency 


($) 


rh 






W(x < x5) W W(x < X;) 






1 IU00 

(3 ) (%) (3) (3) 
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TABLE 7-3 Contanwed 


TRIANGULAR DISTRIBUTION BETA DISTRIBUTION 





Relative Cumulative Relative Cumulative 
= Frequency Frequency Frequency Frequency 
($) 
v= ft W(x as 4) 
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Figure 7-5 
Triangular and Beta Distribution Curves 
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CHARTER ya 
DRY HOLE COSTS 


In Chapter V we discussed the concept of probability 
of success and how it related to the number of successful wells 
to be expected in a drilling program. To account for the num- 
ber of wells that are unsuccessful we have included dry hole, 
Or unsuccessful well, costs into our profitability model. 

While it can be argued that dry hole costs should have 
a continuous frequency distribution, just as the other parameters 
in the model, we will assume for simplicity sake that the dry 
hole costs are constant in this thesis. 

To be consistent with the values developed in Table 
7-1, we will assume that the drilling venture is to take place 
in the Southern Louisiana on-shore area. From January to Dec- 
ember 1970 there were a total of 1245 wells drilled in this 
Mrea witha total footage of 12,703,309 [51]. Therefore, the 
average well depth would be equal to 10,200 feet/well. Using 
1969 average dry hole costs of $13.23/ft [50], the average dry 
hole cost per well would be $135,000. Allowing for inflation 
and deeper wells we will assume a constant dry hole cost of 


$150,000 for each unsuccessful well in the program. 
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CHAPTER IX 


COMPUTER PROGRAM 


The FORTRAN IV computer program, as contained herein, 
was developed to apply the Monte Carlo simulation technique to 
an exploratory well drilling venture using the concepts and 
relationships established in previous chapters, to determine 
the best estimate of discounted profit and to accumulate data 
for each iteration so that a cumulative probability curve can 
be plotted. 

The program'is designed to call up the same series of 
random numbers with each separate computer run so that sen- 
Sitivity studies can be made. The random number generator 
used in this program is a function designed for the CDC 6600 
computer system and is called by RANF(O). The generator pro- 
duces pseudo random numbers, distributed uniformly in the in- 
merval O<x<l. 

The program accepts as data: the number of wells 
drilled, probability of success, dry hole costs, and the value 
M, which is the number of end-points entered in Line 5 plus 
one, (Line 4); the median value of Q, the value of Q at the 
median + 1 standard deviation, Q84, and the boundary conditions, 
Qmax. and Qmin., (Line 6); and the lowest, most likely and 
highest profit per barrel, A, B and C, (Line 7). In addition, 
the end points of the random number ranges corresponding to the 
probability of the number of successful wells is entered as NEP, 
Meine 5). 

(7 
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With this data entered, values are initialized (Lines 
Meelis, 14, 15, 18, 20) and the Standard Deviationeand. Meamecet 
the Lognormal distribution are calculated (Lines 16 and 17). 

For each iteration specified (Line 23) the first series 
of operations establishes the number of successful wells (Lines 
24-32). To do this a random number is called, RANF(0O), mult- 
iplied by 10,000 and compared with the end points corresponding 
to the number of successes’ previously entered (Line 5). Thus, 
the random number called esrabitehoe the number of successes, 
NSW. 

If the number of successes is greater than 0, the pro- 
gram is designed to calculate the size reserve for each success 
and to total them to obtain the value QSUM, which is the total 
reserve size discovered (Lines 37-49). For a detailed descrip- 
tion of the Lognormal distribution generator see Appendix A. 

With the total reserve size determined the profit per 
barrel, for this iteration, E, assuming a triangular distribu- 
tion, is calculated (Lines 51-57). For a detailed discription 
Of the Triangular distribution ganerator see Appendix B. 

Having calculated or inputed all the parameters the 
iteration profit, PROFIT, is computed (Line 58) using the re- 
lationship expressed in Chapter III, Equation 4. Lines 59 and 
60 accumulate this iteration profit and average the value, PTOT, 
over the number of iterations run, MM. XSQ and XSQTOT in Lines 
61 and 62 are used to calculate the standard deviation and coe- 


fficient of variation later in the program. 
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Lines 64 through 176 are inserted to store the values 
Sf iteration profit in arrays by PROFIT intervals in order to 
compute the probability of the interval, PROB, and the cumu- 
lative probability, CUMP (Lines 189-202). 

The variance, VAR, standard deviation, SD, and coeff- 
icient of variation, CV, are calculated for the particular set 
of model parameters in Lines 178, 179 and 180 using information 
previously determined in Lines 61 and 62. 

A listing of the FORTRAN IV computer program is con- 


tained on the following pages. 
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300 


10 
11 


32 


40 


46 
48 


49 


47 
45 


50 


S6 


PROGRAM PROFIT (INPUT s,OUTPUT) 
DIMENSION NEP (25) IY(40)9 QLOGN(20)+ D(40)9 YY(40)s PROB(40) » 


1 CUMP (40) 


READ 100% NePeCUem 

READ 110% (NEP(K)» K=1] 9M) 

READ lll,» QMEDIANs Q84_5 QMAx, QMIN 

READ 1129 As Boe C 

PRINT 15S 

PRINT 113 

PRINT 114+ Neo Ps CUs QMEDIANs Q845 QMAXs QMIN 
PRINT 160% As Bs C 


PkKOF = 0.0 
PROFIT = 0-0 
XSQ = 0.0 


xSQTOT = 0-0 

SO = ALOG(Q84) = ALOG(QMEDIAN) 
TMU = AtLOG(QMEDTAN) 
MM = ] 

PRINT 115 

DO 300 I = 1240 
IY(I) = 0.0 
CONTINUE 

DO 99 II = 141000 
CALL RANF (0) 

Ww] = RANF (0) 

w = wl1#10000. 

J = W 

DO 10 I = leo™ 

IF (JeleeNEP(I)) GO TO 11 
CONTINUE 

NSW = l=) 

IF (NSW.FQ-0-0) 32940 

QSUM = 0.0 

E = 0-0 

GO TO S6 

CONT INUE 

QSUM = 0.0 

00 45 I = leNSW 

RA = RANF (0) 

RB = RANF (0) 

V = (-2.0#AL0G(RA) ) ##e5#C0S (2.083.214) 6#RB) 
QNORM = V#SD «+ TMU 

QLOGN(I) = EXP (QNORM) 

IF (QLOGN(I)-QMIN) 46947948 
QLOGN(I) = QMIN 

GO TO 47 

IF (QLOGN(I)-QMAX) 47947949 
QLOGN(I) = QMAX 

QSUM = QSUM + QLOGN(I) 

CONTINUE 

CALL RANF (0) 

WXLTXI = RANF (0) 

WXGTXIT = 1-0 = WXLTXI 

XI = A # SQRT(WXLTXI* (B-A) *# (CA) ) 
IF (XT.LE.8) 55950 

XI = C = SQRT (WXGTXI#(C-B) #(C-A) ) 
E = X1/1006¢ 

PROFIT = QSUM#E = (N=NSW)#CU 


OON OUS WN = 
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PROF = PROFIT + PROF 


PTOT = PROF/MM 
xSQ = (PROFIT) ##2 


KSQTOT = XSQTOT + XSQ 
QSUM, 


PRINT ITI6« MMe NSWe 
MM = MM ¢« J 

IF (PROFIT. LT e-10-0) 
PeGPROFITsLT. 0-00) 
IF (PROFIT.LT. 1020) 
IF (PROF IT.LT. 20.0) 
IF (PROFIT LT. 30-0) 
IF (PROFIT LT. 400) 
IF (PROFIT.LT. 5000) 
IF (PROFIT.LT. 6040) 
IF (PROFIT LE. 70-0) 
IF (PROFIT LT« 8040) 
IF (PROF IT.LT« 9040) 
IF (PROFIT. .LT.I100.0) 
IF (PROF IT.LT.110.0) 
IF (PROF IT.LI.~120-0) 
paeresOF TT .LI.130.0) 
IF (PROFIT. .LT.140¢0) 
IF (PROF IT.LT-150-0) 
IF (PROF IT.LT160-0) 
IF (PROF IT.LT.170.0) 
IF (PROF IT.LTe180-0) 
IF (PROFIT LI «1906-0) 
IF (PROFIT .LT 200-0) 
IF (PROF IT .LT-210-0) 
IF (PROF ITLL Te 220-0) 
IF (PROF IT.LT 23060) 
IF (PROF IT.LT 240-0) 
IF (PROFIT.LT 250-0) 
IF (PROF IT.LT.260-0) 
IF (PROF IT.LTe270-0) 
IF (PROFIT LI e280-0) 
IF (PROF IT.LF e290¢0) 
IF (PROF IT.LT 300-0) 
If (PROF IT.LT.350-0) 
IF (PROF IT.LT 400-0) 
IF (PROF IT.LT 2450-0) 
IF (PROF IT.LT 500-0) 
IF (PROF IT.LT.S50-0) 
GU TO 99 

Pret? = I1¥(1) + JI 
GO TO 99 

Ty¥(2) = 1Y¥¢2) + I 
GO TO 99 

Ieee = J¥(3) * I 
GO TO 99 

Ty¥(4) = 1y¥(4) + I 
GC 10 99 

Ty¥(5) = 1y¥(5) ¢+ J 
GO TO 99 

1¥(6) = I1Y¥(6) ¢ I 
GO TO 99 

I¥(7) = 1¥(7) + | 
GO T0 99 


GO 
GO 
GO 
GO 
GO 
GO 
GO 
GO 
GO 
GO 
GO 
GO 
GC 
GC 
(,0 


GO 


TO 
TO 
TO 
TO 
TO 
TO 
TO 
TO 
TO 
TO 
TO 
TO 
TO 
TO 
TO 
TO 
TU 
TO 
TO 
TO 
TO 
TO 
10 
TO 
TO 
TO 
TO 
TO 
TO 
TO 
TU 
TO 
TO 
TO 
TO 
10 
TO 
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PROF ITs 
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79 
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Bg 
90 
91 
G2 
93 
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96 


IY(8) = 
oO FO 9S 
Iy(S) = 
GU TO 99 
IY¥(10) = 
GO TO 99 
ly(ll) = 
GC 10 99 
ly(le) = 
GO 10 9S 
Prtl3) = 
mG TO GJ 
[y(ls#) = 
GO TO 99 
metcisS) = 
ov TO 9S 
Prtil6) = 
WO TO 99 
mats) = 
wl FTO ¥9J 
my tls) = 
nu TO 99 
PYti9S) = 
GO TO 99 
Iy(ed) = 
GC TO 99 
Hricly = 
“Cc 10 9S 
l¥(2e) = 
130 TO 99 
ly(23) = 
GC 10 99 
Iy(24) = 
GO TO ¥9g 
Iy¥(2eS) = 
GC TO 99S 
ly(26) = 
UU TO 99 
Pyae?) = 
GU TU 99 
ly(28) = 
GU TO 9S 
ly(?9) = 
GO TQ 99 
tY¥(30) = 
GO TO 99 
TY¥(31) = 
GQ TO 99 
Iy(32) = 
GO TO 99 
Iy¥(33) = 
GC TO 99 
IY¥(34) = 
Go TO 99 
Iy(35) = 
GC TO 99 
Ily(36) = 
GO TO 99 


TY (8) 


TY (9) 


1Y¥ (10) 
TY (11) 
1Y( 12) 
1Y¥(13) 
TY (14) 
1y¥ (15) 
Ty (16) 
LY ¢rh) 
TY (18) 
1y¥ (19) 
1Y¥ (20) 
LY 12) 
ete? 
lY (23) 
TY (24) 
IY (25) 
1Y (26) 
rte?) 
ly (2e) 
LY (29) 
TY (30) 
TY (31) 
TY (32) 
1Y¥(33) 
1Y (34) 
ly (35) 


LY (36) 
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220 
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200 
100 
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140 
141 


144 
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160 


lY¥(37) = I1Y¥(37) ¢« ] 

CONT INUE 

MM = MM = ] 

VAR = (XSQTOT = (MM#(PTOT##2)))/(MM = 1) 
SD = SURT (VAR) 

CV = SD/PTOT 

PKINT 140° MMoPTOT 

PRINT 141% VARs SDoe CV 

PRINT 144 

PRINT 155 

PRINT 113 

PRINT 1149 Ns Ps CUs QMEDIANs Q845 QMAXs QMIN 
PRINT 160% As Bs C 


PRINT 145 
SUMM = 0.0 
el). = -10.0 


DO 200 IM = 1.437 

YY(IM) = IYCIM) 

PROS (TIM) = YY(IM)/1000- 

SUMM = SUMM ¢ PRO WQ(IM) 

CUMP (IM) = SUMM 

Ky = 1M + 1 

1F(IM-32) 210+ 220+ 220 

D(KJ) = DOC(IM) + 10-0 

GO TO 230 

D(KJ) = O( IM) + 50.20 

CONT INUE 

PKINT 150+ OC(IM) + PROR( JIM) » CUMP(IM) 

CONT TINUE 

FORMAT (10X91S92F1023e15) 

FORMAL (10X%01714) 

FOKMAT (10X04F10.3) 

FORMAT (10X93F 10-2) 

FCRMAT (///56X9*NOe OF WELLS*93X%5#PROB OF #9 2X e9#DRY KOLE%9/99Xe#ORI 
ILLED* 95X 9 *#SUCCESS* 9 4X 9 FCOST #9 4X9 #QMEDI AN* 9 3X 9 #084 %95X 9 FQMAX* OSX 
2*QOMIN*) 

FORMAT (1OXe 1S 8X oF Ge394XaFGe3d 9 IXoF Te 2alXoF Tele 2XoF lelalXaF le2//) 

FORMAT (Sx5*# I TERATION®, 3X 9#NUMBER OF *94Xo9#F IELD#,2X9#PROF IT PER#,2Xx 
Le®ITERATION*+3Xe*#BEST ESTIMATE*# 9/1 7X9*SUCCESSES# 94Xe*#SIZES*9SXo*#BA 
CRREL * 96X99 #PROF I] T*e+6Xe#0F PROF IT*# s//) 

FORMAT ( SXeIlGsll1X913 s6XeF Gel sbXaF4e2eXoF Gel a BX oF 602) 

FORMAT (///e5Xe*THE BEST ESTIMATE OF PROFIT BY THE MONTE CARLO SIM 
LULATION METHOD*#eo/7Xe*#AT %e9I14(% ITERATIONS IS oF 7.29% MILLION DOLL 
1AKS#*#) 

FCRMAT (/// 95X98 THE VARIANCE IS oF 12¢39/95Xe* THE STANDARD DEVIATI 
1ON [S *oF] 2e639/9SXe*THE COEFFICIENT OF VARIATION IS *9F 6.4) 

FORMAT (1H1) 

FORMAT (15Xe#UPPER LIMIT OF #9 7Xe#PROBABILITY#+5X9*CUMMULATIVE*#, 
17e1S5Xe*#PROFIT INTERVAL *#.S5Xe2#0F INTERVAL*®# 9SX 9 #PROBABILITY*—9///) 

FORMAT (20X oF Ge0915X FS e491 1X 9F 504) 

FORMAT (lH1]91SXe#PROFIT DISTRIBUTION OF AN EXPLORATORY ORILLING PR 
1OGRAM#) 

FORMAT (/15X-e#THE TRIANGULAR OISTRIBUTION FOR DETERMINATION OF THE 
l*9/15Xe#PROFIT PER BARREL HAS THE FOLLOWING EXTREMES IN CENTS <— #, 
C7 9lLSXeFA = *oFGe2oTX9*#B = Ref Bele TXo*#C = *oF6e2s///) 

END 
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CHAPTER x 


BEol ESTIMATE OF DISCOUNTED EXPECTED PROFIT UstNGepAoeee 


In previous chapters the profitability model was dev- 
eloped and the individual parameters and their frequency dis- 
tributions were discussed. In addition, certain Base Data was 
established for each parameter so that the best estimate of 
saescounted expected profit, E(P), could be calculated. This 


Base Data 1s summarized in Table 10-1. 


TABLE 10-1 
BASE DATA 
Parameter Distribution Value 
Number of Wells Drilled - n : 20 
Probability of Success - p Binomial Biles 
Reservoir Size Lognormal 
Q at 5% probability 5 x 10° bbl 
Q at 95% probability 79 X i0° bbl 
Q maximum 140 x 10° bbl 
Q minimum 1 x 10° bbl 
Discounted Profit/Barrel dete eablielie 
Lowest - e, $. «5.0/ bbil 
Most bikely - e, $ .88/bb1 
bagihest - e€. $1.15/bbl 
Dry Hole Cost - C,, None IL S05 OU, 
Number of Iterations ESOog 


8 4 
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On the basis of this Base Data and using the Fortran 
IV program included in Chapter IX a value of Best Estimate of 
Discounted Expected Profit, E(P), of $68.02 million was obtained. 
See Appendix D for a printout of the number of successes, 
fecerve Sizes, profit per barrel, lteration profit and average 
best estimate of discounted profit for each of the 1000 iter- 
ations using the Base Data. In addition, Appendix D contains 
cumulative probability data for each specified interval of E(P). 
Using this cumulative probability information the histogram 
included as Figure 10-1 was prepared from which it can be seen 
that: 
1. The Median Value (at 50% probability) of Discounted 
Profit is $56.03 million. 
2. With 90% confidence we could expect a discounted 
profit of between $1.51 million and $170.0 million. 
3. The probability of a loss is 4.2%, or stating it 


another way, there is a 95.8% chance of a profit. 


Further, from the printout in Appendix D we see that: 
4. The standard deviation is $53.66 million. 
5. The coefficient of variation is ./7889. 
The significance of the standard deviation and coefficient of 


variation as measures of risk were discussed in Chapter IV. 
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Figure 10-1 
Cumulative Probability - Base Data 
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CHAPTER Gl 


SENSITIVITY ANALYSIS 


Freedom from uncertainty is a luxury rarely enjoyed by 
the contemporary decision maker. In a society characterized by 
change, uncertainty is an established and accepted fact of life. 
This is especially true in the oil business. Notwithstanding 
the difficulties, the decision maker is expected to attain 
satisfactory results as measured by the goals of his organiza- 
ion. 

In the face of uncertainty, the most recurring quest- 
ions he faces are of the form, "What if...? 

1. What if labor costs rise and dry holes cost twice 

as much? Three times as much? 

2. What if the posted price of oil decreases by 30%? 

3. What if we drill 30 exploratory wells rather than 

20? How about 10 wells? 
The "what if'' question may be viewed as an introduction to ths 
sensitivity analysis. In its mathematical sense, sensitivity 
analysis is a study to determine how possible changes or errors 
in parameter values affect model outputs. In its applied or- 
ganizational setting, sensitivity analysis may be broadly de- 
fined as a study to determine the responsiveness of the con- 
clusions of an analysis to changes or errors in parameter 
values used in the analysis. In other words, would our decision 
be any different if a parameter value, such as the probability 


Be success insout profitability model, were to change. The 


87 





ee OO 


a 





88 


prime function of sensitivity analysis can be said to be to 
meovide a better understanding of risk. Sensitivity analysis 
eam provide guidelines for allocating organizational resources 
to data collection and data refinement activities [40]. The 
benefits of subjecting models to sensitivity analysis are 
summarized by Arnoff and Netzorg [2]: 

The use of an operations research model is especially 

important and advantageous in that it permits experi- 

mentation "on paper", without manipulation of the act- 
ual system. In using the model, one can assess the 
sensitivity (response) of the system to a wide variety 
of conditions - without requiring either the time, 
expense, or risks associated with experimenting with 
the system itself (if such experimentation would, in 
fact, be possible and meaningful). Hidden relation- 

Ships can be brought to light and brought to bear upon 

ScocisiOnsS and Contnole Of aetivity, 

In this thesis a profitability model for an exploratory 
well drilling venture was developed and tested with Base para- 
meters. A decision problem is now faced. If the model was 
developed in order to provide results which would serve as a 
basis for recommending action to management, it is possible, 
me this point, te use the model results as 1s, or to select 
miner alternatives: to collect additional data: refine initial 
parameters: or even alter the model itself. On the significance 
Mee this decision, Willian T. Morris [35] comments: 

The nature of management science is such that sooner or 

later it has to recommend some action other than making 

more observations. Usually this comes sooner rather 


than later. Eventually it will suggest that action be 
taken consistent with either acceptence or rejection 


of the hypothesis. It will recommend management "act 
as if" the hypothesis were true or "as if" it were 
false. Indeed, the fundamental question in management 


science may well be when to stop collecting data and 
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developing the model, and when to produce a recommen- 
dat Lon = ore aiG Eaton. 


In the past the emphasis has been on whether information re- 
quired for decision making could be developed. With our grow- 
ing capacity to generate information or data, the emphasis now 
appears to be shifting to the question as to what extent this 
Capacity can be utilized beneficially. It may be recognized 
Haat while the collection and refinement of data are costly, 
it may serve to improve recommentations for action. In most 
cases, however, little basis exists for projecting whether the 
additional information is worth the cost. Statistical decision 
theory, and a model such as the one developed in this thesis, 
will assist in the decision as to whether additional informa- 
Eion is "worth the cost". The basic requirements for the stat- 
istical decision theory model are summarized by Morris [36] as 
follows: 
[The Model] must be able to enumerate the possible 
outcomes of future data collection efforts and, further, 
moncompute the probabilities of these outcomes. In 
addition, it must be possible to indicate just how the 
information willquantitatively change the decision 
maker's view of his choice. 
Mmiean excellent example, which illustrates this concept, C. 
Jackson Grayson [21], demonstrates the application of decision 
theory to drilling decisions by oil and gas operators and also 
to the decision problem of whether to purchase additional in- 
formation before making a final drilling decision. 


Clearly then, a sensitiviey analysis which discloses 


decision insensitivity to estimated parameter variations 
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mmeecludes the necessity Of@antroducing more costly. time 
consuming data, collection data refinement and statistical 
decision models into the analysis. 

It is in this vein that the sensitiviey studies are 
conducted in this thesis. The profitability model provides 
Bieemecessary vehicle for the analysis. From the results it 
1s possible to determine the outcomes, as well as probabilities 
of these outcomes, associated with changes in parameter values 
thus providing some answers to the "What if?'' questions posed 


moeene beginning of this chapter. 





CHAPTER XII 


SENSITIVITY STUDaaes 


ine Outline of Studies 

The sensitivity studies were conducted using the 
profitability model developed in Chapter III. In these studies 
each parameter in the model was varied from the base values as 
summarized in Table 10-1. Each time a parameter was varied all 
remaining parameters were set to their base values and the FOR- 
TRAN IV program, developed in Chapter IX, was executed 1000 
times. In this manner a printout similar to the one for the 
Base data, included in Appendix D, was obtained for each para- 
meter change or Case. 

Thus, for each Case, at 1000 iterations, we have the 
best estimate of discounted profit, the standard deviation, the 
coefficient of variation, the iteration profit, cumulative 
probability data, etc. Tables 12-1 through 12-9, grouped by 
parameter being varied, contain a summary of some of this data 
(at 1000 iterations) as well as extrapolated information such 
as the median discounted profit, the maximum and minimum dis- 
counted profits within 90% confidence limits and the chance 
of making a profit greater than zero. The values using the 
Base data have been included on each Table to assist in eval- 
uation. It should be noted that each time the model was run 
the same series of random numbers were called up so that com- 
parisons between Cases could be made. Figures 12-1 through 
12-8 contain histograms for several Cases, plotting the cum- 


sill 
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ulative probability for each Specified interval of discounted 
expected profit. Figures 12-1 through 12-8 are similar in 
presentation to Figure 10-1 for the Base data. On each hist- 
Oogram the 90% Confidence Limits have been represented by dashed 
lines at 5% and 95% cumulative Probable tiess 

To assist in the comparison of each Case to the Base, 
mepres 12-1A through 12-9A have been prepared and are inserted 
behind the Table it serves to amplify. These Tables contain 
the percentage differences between the Case parameter set 
values and the Base parameter set values. The percentages 


were calculated using the following relationship: 


(Case data value) - (Base data value) 


Base data value x 100, (1) 


A(in %) = 


In addition to the Cases noted above, the Base data 
was run for 10,000 iterations rather than 1000 iterations. 
Additionally, the Beta Distribution Variable Generator was 
substituted for the triangular distribution of profit/barrel 
in the FORTRAN IV Program and the Base data Was rune tor 000 
iterations. The results of both these Variations are summar- 


ized on Table 12-10. 


12.2 General Observations 
pee Se Vations 
One of the prime purposes of the sensitivity studies 
is to assess the change in observed values meSULtIng frome 


Change in parameter value. The majority of the effects of 
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Figure 12-3 


Cumulative Probab 
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1lity - Q at 95% Parameter Set 


Figure 12-4 
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Figure 12-5 
Cumulative Probability - Ch Parameter Set 
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Figure 12-6 
Cumulative Probability - e(a) Parameter Set 
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Figure 12-7 
Cumulative Probability - e(b) Parameter Set 
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these parameter changes can be observed on Tables 12-1 through 
12-9 and 12-1A through 12-9A and require no further datailed 
discussion. The analysis of a parameter change to a change in 
Best Estimate of Discounted Profit and to the risk as measured 
by the changes in Coefficient of Variation will be discussed 
im the next Section. 

At this point however, some general observations 
regarding the sensitivity of the model to parameter changes 
can be made. First, an inspection or the Tables will reveal 
that the one single observed value which fluctuated the great- 
est on a percentage basis was the minimum value of discounted 
profit which could be expected with 90% confidence. The reason 
for this high percentage fluctuation was the relatively small 
numbers involved. Thus, a variation of $7.87 million on a 
$1.51 million Base resulted in a percentage change of approx- 
imately 980%. 

The second observetion made from the Tables is that, 
except in isolated cases, the chance of a profit greater than 
zero was fairly constant at approximately 95 to 96%. Stating 
it another way, there appears to be only a 4 to 5% chance of a 
loss in a majority of the Cases. When the probability of 
success was reduced to 1 in 20 there was still a 62% chance of 
making a profit and, even when the number of wells drilled was 
reduced to 5 rather than 20, the probability of a profit was 
approximately 56% (the lowest value in any Case). 


An inspection of Table 12-10 reveals that themsubstieu- 
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tion of the Beta distribution for discounted profit per barrel 
rather than the Triangular distribution has very little effect 
on the Best Estimate of Discounted Profit, E(P). The risk as 
measured by the Standard Deviation and the Coefficient of 
Variation is increased however, by about 4.4%. The minimum 
value of discounted profit with 90% Confidence Limits has been 
decreased by 4.34% and the maximum value has been increased by 
1.78%. Additionally, the range of values with 90% Confidence 
Limits has been increased with the Beta distribution by 2.24%. 
With regard to the Base data run at 10,000 iterations 
rather than 1000 iterations, it is seen from Table 12-10 that 
the value of E(P) only varies by 1.47% from the 1000 iteration 
value while the relative risk as measured by the Coefficient 
of Variation increases by 1.5% primarily due to the decrease 


in the absolute value of E(P). 


teow artects of Parameter Changes on the Best Estimate of 


Discounted Profit and Coefficient of Variation 

The first question explored in this Section is the 
effect of a change in parameter value on the best estimate of 
discounted profit, E(P). This effect shall be assessed by 
calculating a Significance Ratio for each Case. The Signifi- 
cance Ratio is defined as the change in observed value (i.e. 
E(P)) to a change in the parameter value. In a conceptual 
sense the Significance Ratio measures the significance of a 
change in a parameter of the model on a value calculated by 


the model. We assume, of course, that all other parameters 


ZZ 


(except the one being studied) remain the same. Mathematically, 


this can be expressed as: 


A Observed Value 


SR = XK Parameter (1) 


With this Significance Ratio the relative importance of a 
change in various model parameters on the value being observed 
can be measured. For the profitability model it is possible, 
for example, to measure the relative importance of a change 
mcr y hole costs, lowest profit per barrel, etc. on the best 
estimate of discounted profit for the total venture. For the 
analysis of the best estimate of discounted profit, Equation 


1 can be rewritten: 


ete) 
SR = A Parameter (2) 


The values to be inserted in Equation 2 are obtained from 
Tables 12-1A through 12-9A. As an example, from Table 12-1A 
the following Significance Ratio is obtained: 

For a change in parameter of -75.00% (Case 1) from the 
Base we have a commenserate reduction of E(P) of -76.77%. 


Inserting these values in Equation 2 we have: 


~ ne a5 
SR = | Sys 00) | 
Sh 3.02 


For Case 2 we have the following relationship: 


-50.44 
ee 


= 05,00 
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In a similar fashion the Significance Ratios can be calculated 
for each of the nine Cases involving a variation of the number 
@z wells drilled as shown on Table 12-1 and 12-14. (lf ethics 
done Significance Ratios ranging from .81 to 1.18 with a mean 
value of 0.99 are obtained. 

Table 12-11 contains a summary of the Ranges and Mean 
Significance Ratios for each of the parameter sets contained on 
Tables 12-1 through 12-9 and Tables 12-1A through 12-9A. In 
addition, Table 12-11 contains a value of the Normalized 
Significance Ratio which reduces each Significance Ratio to a 
value on a 1-100 scale in proportion to its actual Significance 
Ratio. The Normalized Significance Ratio can be utilized to 
facilitate comparison of parameter sets 

TABIGE 2 = 
SIGNIFICANCE RATIOS - E(P) 


Significance Ratios 


Parameter Set Range Mean Normalized Mean 

n ol ewes .99 96 

Pp #916: -  lyeko Leads 100 
Qat 5% .19 - .40 £50 29 
Q at 95% SSomas..82 ae 70 
Q (max) NO 2 -= eee sOLT OF 
Q (min) 0 0 0 
os 0 .04 04 
e (a) 0 peu 20 
e(b) 0 . 36 GS 


ec} 0 47 46 
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several observations can be made from Table 12-11. 
First, it can be seen that a change in the number of wells 
drilled and probability of success are both almost equally 
Significant as far as their effects on the value of E(P). In 
fact, there is virtually a 1 to 1 relationship between a change 
in the number of wells drilled and probability of success and a 
change in the best estimate of discounted profit. In other 
words, if the number of wells drilled is increased by 10%, 
almost a 10% increase in discounted profit can be expected. 
Likewise, if the probability of success is decreased by 20% 
about a 20% decrease in the discounted profit can be expected. 

Also seen from Table 12-11 is that the next most 
significant parameter in the model is the estimated value of 
reserve size with 95% probability, Q at 95%. In fact, this 
parameter has a normalized significance of 70 which is approx- 
imately 2.5 times as critical as the lower estimate of the 
Wememve size, Q at 5% probability, on the discounted profit. 
In a similar fashion it can be observed that variations in the 
values of the upper and lower boundary limits for reserve sizes 
faemlittle or no significance on the discounted prorit. 

Although one might intuitively feel that changes in dry 
hole costs would have a material effect on the discounted 
profit it is seen that it only has a normalized Signimreance 
of 4 which, in essence, is relatively insignificant in compar- 
ison to the other parameters. 


In comparing the three parameters associated with the 


discounted profit per barrel it is noted that, of the three, 
changes in the highest profit per barrel, e(c), are 2.3 times 
aS Significant as changes in the lowest value, e(a), and 1.3 
times as significant as changes in the most likely value, e(b). 
Also observed is that of the 10 parameters measured, the value 


@imeyc) ranks fourth in significance. 


A similar analysis regarding risk can be made by 
computing the Significance Ratios for the Coefficient of 


Variation. In this case Equation 1 can be rewritten as: 


Aue 
aS A parameter (3) 


Table 12-12 contains the results of these calculations. 
TABLE 12-12 
SIGNIFICANCE RATIOS - CV 


Significance Ratios 


Parameter Set Range Mean Normalized Mean 
n Gish IES Ie 167 100 
p allie) , So 64 96 
Q at 5% ,046> 43 qo. 28 
Q at 95% /01 - .18 . 06 9 
Q (max) mOOMaee: 15 .10 ies 
Q(min) 0 0 0 
Cc, pds 05 . 04 6 
e(a) BOG = wecler .09 ile: 
e(b) 0 - .06 nS 5 
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From Table 12-12 certain observations regarding the 
Significance of a change in model parameters on the relative 
risk of the venture can be made. First, it is observed that a 
change in the number of wells drilled and in the probability 
of success both have about the same effect on the risk of the 
venture. It is interesting to note that the Normalized 
Significance Ratios for E(P) are exactly reversed for CV. In 
this regard, Table 12-11 indicated that a change in Q at 95% 
had a greater significance than a change in Q at 5%. Table 
1Z-1Z reveals the opposite is true for the relative risk. 

Here a change in Q at 5% has a much greater significance than 
a change in Q at 95%. The same reversal is evident in the 
values of lowest and highest profit per barrel, e(a) and e(c). 

It must be remembered, of course, that the Signifi- 
cance Ratio for risk (coefficient of Variation) does not in- 
dicate whether the risk is increasing or decreasing with a 
parameter change, but only the relative amounts that it is 
Changing. For the increase or decrease correlation we must 


g0uto Tables 12-1 through 12-10 and 12-1A through 12-9A. 


eat tonvergsencesot the Expected Profit 


If the Monte Carlo simulation model is executed only 
once, only one value for the discounted profit is generated and 
there is no distribution of values that can be analyzed. Since 
Simulation is being used because the model is too complex for 
other types of analyses, there is no basis for relying on the 


single result as a typical output of the model - there is no 





Ub, 7 


way to tell whether the result is an extreme case or about the 
average. 

As the simulation process is repeated frequencies of 
occurrence are obtained from which the averages, the extremes, 
and what the shape of the distribution would be if the model 
had been executed infinitely many times can be inferred. In 
this regard, R.F. Barton [6] poses the following question: 

how reliable are the results from a limited 

number of executions? Two executions may give us a 

great deal more information than a single execution. 

One hundred iterations may give us hints of the extre- 

mes and shapes of the ultimate distributions. A thou- 

sand executions may perhaps give us results that app- 

roximate smooth curves. But a second thousand execu- 

tions may suggest smooth curves of different locations 
and shapes. A third thousand executions may generate 

results that look even different. On the other hand, 

the separate sets of one thousand executions may look 

so similar we feel no need for further executions. 

These possibilities raise a very practical question: 
How long to run a simulation model? One aspect of this quest- 
ion was analyzed by measuring the convergence of the discounted 
Expected Profit,E(P), averaged over the number of iterations 
run, with the discounted Expected Profit averaged over 1000 
iterations. 

Ncean example, if “the valluewEQP)Meer the famst Peeration 
is $100¥0 million and $50.0 million for the second iteration 
the discounted E(P) as determined by two iterations is $75.0 
million. We can then compare this value with the E(P) averaged 


over 1000 iterations, say $60.0 million, to determine the 


rapidity of convergence. If we were satisfied with an Expected 
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Profit that varied no more than + 5% from the value which would 
be obtained with 1000 iterations (90% certainty) then by plot- 
ting or analyzing the data we could determine how many iter- 
ations were necessary to reach this point. It may turn oyblic aly) 
our example, that a value of $60.0 million + $3.0 million was 
achieved after only 500 iterations. By performing this analysis 
with many different values of model variables, using the same 
random numbers, and each time determining the point beyond 
which the averaged E(P) varies by no more than + 5% of the 1000 
iteration value, we may be able to infer the number of itera- 
tions necessary for 90% confidence limits for this particular 
model. 

This is, in fact, what has been done. For each set of 
model variables the iteration beyond which the value of E(P) 
does not vary by more than + 5% from the value of discounted 
expected profit after 1000 iterations was measured. These 
Values are contained on Table 12-13. Figures 12-9 through 
12-14 graphically depict this convergence process for various 
combinations of model parameters. 

In addition to this analysis, the model was run with 
the base data for a total of 10,000 continuous iterations to 
determine if any periodicities in the variables, as suggested 
by Barton [7], occur which would materially effect the value 
of the discounted Expected Profit after 10,000 iterations. 
Tate 12-10 and Figure 12-15 contain the results of this test- 


It can be seen from Figure 12-15 that no periodicities occur 
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and the 10,000 iteration value of E(P) is, in fact, only 
slightly different from the 1000 iteration value. 

A visual inspection of the data contained on Table 
12-13 indicated that Expected Profit, E(P), does not converge 
with any regularity, either within 80% or 90% confidence limits. 
It is seen that even within a parameter set(i.e. number of 
wells drilled) the number of iterations vary greatly. The only 
Significant fact which can be discerned from the data is that 
values of E(P) within 80% confidence limits can be achieved 
anywhere from 10 to 339 iterations and within 90% confidence 
limits anywhere from 63 to 710 iterations. If the highest 
and lowest iteration value is eliminated the range for 80% 
confidence limits is 18 to 242 iterations and the range for 


90% confidence limits is 68 to 596 iterations. 
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12.5 Changes in Return and Risk with the Number of Wells 
Drilled 


Since the risk of the venture has been measured by the 
coefficient of the variation it may be of interest to ascertain 
whether there is a mathematical relationship between the risk 
and the number of wells drilled. To determine whether a mathe- 


matical relationship exists we start with Equation 7, Chapter V: 


C= WD 
P 


Se 


Now, if p is held constant it can be seen that the Coefficient 
of Variation should be proportional to the number of wells 


drilled as follows: 


~ gal 4 
CV a= (4) 
ore 


cV x Yn = Constant ‘e>) 


Table 12-14 contains a test of this proportionality hypothesis 


Using the data contained on Table 12-1 
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TABLE 12-14 
PROPORTIONALITY TEST - CV to n 


n Yn Coeft: of Vananion Constant 

5 RITES 1.59 oo 
eG) ble Hall S208 
eZ 3.46 Oe 50 
iS Ste i OR oz 54150 
Ihee 4.24 0.84 Soe 
20 4.47 0.79 Sars 
ne 4.69 ORS Sac 
DO 00 OF Wk Sa 
30 5.48 0.64 54 30 
J oe 92 OAGZ 5.07 


It can be seen from Table 12-14 that the Coefficient of Vari- 
ation is indeed proportional to 1/Yn. Thus, if the number of 
wells is increased by 4 times the risk will be reduced by a 


mieeer of 2. 


It is also interesting to note that the Best Estimate 
of Discounted Profit per well, or E(P)/n, remains fairly con- 
Stant while the risk, as measured by the Coefficient of Vari- 
ation, decreases with 1//n. This can be seen on Table 12-15. 
This would tend to indicate that while the best estimate of 
return per well remains fairly constant with the number of 


wells drilled the risk associated with that return deminishes 





as we increase the number drilled. 
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CHAPTER XIII 


CONCLUSIONS 


fet General 

The following general conclusions are drawn from the 
mesults of this study: 

1. The model developed herein provides a suitable 
vehicle for measuring the sensitivity of the parameters affect- 
ing the profitability of an exploratory oil drilling venture. 

2. The parameters in the profitability model should 
include not only single best estimate values but, where app- 
ropriate, also a probability distribution of values. 

3. The probability distribution best representing the 
number of successes is the Binomial distribution, especially 
when the probability of success and number of wells drilled are 
in the order of magnitude normally encountered in an exploratory 
drilling program. For probabilities of success approaching 50% 
and large numbers of wells drilled, the Normal or Gaussian dis- 
tribution can closely approximate the Binomial distribution. 

4. The reserve size discovered by a successful well is 
best represented by a Lognormal distribution. The values can 
be calculated manually using points scaled from a straight line 
drawn on log-probability paper between values of reserve size 
which can be expected with 5 and 95% probability. The lognormal 
distribution can be simulated in computer applications with 
excellent correlation using the Lognormally Distributed Variable 
Generator contained in this study. 


az 





lene 


>. The profit per barrel received over a period of 
years from oil produced must be discounted to reflect the time 
value of money. 

6. The Triangular distribution function adequately 
serves to represent the discounted profit per barrel. The Beta 
distribution function provides only slightly different values 
of discounted profit per barrel when tested using the same 
input data used in the Triangular distribution. Both the 
Triangular and Beta distributions can be simulated in computer 
applications using the Generators contained herein. 

7. The complexity of the distributions in the model 
precludes manual computation of profitability and sensitivity 
studies thereby necessitating the use of the Monte Carlo 


technique and a digital computer. 


MomammocnsSitivity Studies 


In the sensitivity studies, where individual parameters 
in the model were varied, with the remaining parameters being 
held constant, the following conclusions were drawn: 

1. The one single observed value which fluctuated the 
greatest (on a percentage basis) was the minimum value of 
discounted profit within 90% Confidence Limits. 

2. Except in isolated cases, the chance if a profit 
greater than zero in all tests remained fairly constant at 95 
to 96%. Conversely, there remained a fairly constant chance 


of 4 to 5% of incurring a loss. 
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3. <A sSubsitiution of the Beta distribution for dis- 
counted profit per barrel rather than the Triangular distri- 
bution makes very little effect on the best estimate of dis- 
counted profit for the venture or the risk as measured by the 
coefficient of variation. 

4. There appears to bDégno periodicitues mine theapara 
meters of the model and the best estimate of discounted profit 
calculated after 10,000 executions varies only slightly from 
the value calculated after 1,000 iterations. 

5. it is possible to calculate values of the best 
estimate of discounted profit for the venture within 80% 
confidence limits (+ 10% of the 1,000 iteration value of E(P)) 
by executing the model between 18 and 242 times; and within 
90% confidence limits by executing the model between 68 and 
596 iterations. Further, there appears to be no regular 
convergence of E(P) with the number of iterations performed. 

6. <A Significance Ratio, relating a change in Observed 
value calculated by the model to a change in parameter value, 
provides an excellent tool for comparing the effects of changes 
of different model parameters, especially when the Significance 
Ratio is normalized on a 1 to 100 scale. 

7. Using the Significance Ratios to relate parameter 
changes to discounted venture profitability, E(P) 1t- @seeon- 
cluded that with the parameter values used: 

a. Changes in probability of success and number 
of wells drilled have almost equal significance on changes in 


venture profitability. 
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b. Changes in the lower range of reserve size have 
only about 30% of the effect of changes in p orn on E(P). 
However, changes in the upper range have a 70% significance. 
The values selected at lower and upper boundaries for reserve 
sizes in the model have little or no effect on the value or E(P). 

c. Changes in the highest value of discounted 
profit per barrel have more than twice the significance as 
changes in the lowest values of discounted profit per barrel; 
but even then, have only 46% of the significance as changes 
immeeor pon E(P). 

d. Changes in dry hole cost have relatively 
little effect on E(P). 

8. Using the Significance Ratios to relate parameter 
changes to the relative risk of the venture as measured by the 
coefficient of variation, CV, it is concluded that with the 
parameter values used: 

a. Changes in the number of wells drilled, n, and 
the probability of success, p, have about the same effect on 
the risk of the venture. 

b. Changes in the lower range of reserve size are 
next in significance. Changes in this parameter are more than 
3 times as significant as changes in the upper ranges of 
reserve sizes. 

c. The value selected as the upper boundary for 
reserve size will effect the risk in the venture but the lower 


boundary limit will not. 
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d. Changes in dry hole costs do little to change 
the risk of the venture. 

e. Changes in the lowest value of discounted 
profit per barrel is about 1.5 times as significant as changes 
in the highest value of discounted profit per barrel but even 
then, has only a significance of 13% as compared with the risk 
associated with changes in the number of wells drilled. 

9. While the return per well remains fairly constant 
as the number of wells increase, the relative risk of the 
venture decreases as the square root of the number of wells 


drilled. 


The seemingly large number of findings drawn from the 
results of this study are indicative of the value of this type 
Stmanalysis. It is recognized, of course, that not all the 
conclusions listed herein will be considered of equal import- 
ance by the reader - it would be rare if they did. The fact 
that these findings could be drawn at all, however, leads to 
the overall conclusion that sensitivity analysis of this type 


provides a valuable tool for the decision-maker. 





CHAPTER XIV 
RECOMMENDATIONS 


It is in the light of the problems facing most decision- 
makers, namely the "What if...?'' questions cited in Chapter XI 
and the contention so lucidly expressed by William T. Morris 
[35] that at some point we must stop gathering data and re- 
fining our information and start recommending action, that the 
profitability model is ideally suited. 

It is therefore recommended that a model of this type 
be used not to obtain absolute data on whether a profit will 
be made for the venture (or some other observation of this 
type) since other methods of analysis may be better suited for 
this purpose. Instead, it is recommended that the model be 
used to assess the relative sensitivity of the results to para- 
meter changes in the model. With this knowledge of parameter 
sensitivity we can now devote time and effort to refining only 
those parameters which have a significant affect on the outcome. 

For example, if we were concerned with the best estimate 
of discounted profit for the venture we would probably devote 
our efforts to refining our estimates of the probability of 
success, the upper range of reserve size and the highest value 
of discounted profit per barrel (assuming, of course, that the 
number of wells drilled is fixed). We would not, in this 
example, concern ourselves with such items as the dry hole 


costs, the lower range of reserve sizes and the like. 
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In this way, using the sensitivity analysis approach, 
the decision-maker can achieve the maximum return for the 


effort he can afford to expend. 

















THE LOGNORMALLY DISTRIBUTED VARIABLE GENERATOR 


The FORTRAN IV Computer program for lognormally 
distributed deviates requires two different, uniformally 
distributed random fractions to generate one normal variable. 
With two uniformally distributed fractions, RA and RB, one 
generates a variable, V, which is normally distributed about 


a mean of O, with a standard deviation of 1, thus, 


V = (-2 In RA)*? x COS 27 (RB) 


To convert V to a variable QNORM, which is normally distrib- 
uted about a mean of TMU, with a standard deviation SD, we 


Eeiipay use the relationship, 
QNORM = (V) SD + TMU 


By multiplying the natural logarithm base, e, to the normal 
deviate power, a lognormally distributed deviate, QLOGN, is 


generated. 
QLOGN = e QNORM 


Once the Lognormally distributed deviate is generated 
it is tested against the lower and upper limils established as 
inputs to the program. If the deviate is less than the minimum 
value of QMIN, it is set to read QMIN. If it is greater than 
the maximum value of QMAX, it is set to read QMAX. 


It will be noted that the program requires a value of 





mean and standard deviation. These values can be obtained 
from the median value of Q at W( £Q) = 50%, and the value of 


Q at WC SQ) = 84% (or, Y mel Oy) as follows: 


From Equation 3, 
ye (median) = ye (mean) 
From Equation 4B, 
In xX, = We 
Qmedian = Xm 


and, Thu ="y 5 


then, 
TMU = log, (Qnedian) 


TMU = ALOG(QMEDIAN) 


Vemma i 


ey) oy 


ae 


Yn * Sy = In Q at Wg) = 84 
Ol 
yyy tagee= In Q(84) 
and, 
A enn er 
Therefore, 
(Yq * Sy) - Ym = In Q(84) - In(Q,) 
Oe ¢ 


SD = ALOG(Q84) - ALOG(QMEDIAN) 





The FORTRAN IV Computer program incorporating these concepts 


is, listed as follows: 


PROGRAM QLOGN (INPUT ,OUTPUT) 

READ 100, QMEDIAN, Q84, QMAX, QMIN 

SD = ALOG(Q84) - ALOG(QMEDIAN) 

TMU = ALOG(QMEDIAN) 

RA = RANF (0) 

RB = RANF (0) 

V = (-2.0*ALOG(RA)**.5* COS(2.0*3.1416*RB) 
QNORM = V*SD + TMU 


QLOGN = EXP (QNORM) 

IF (QLOGN - QMIN) 6, 7, 8 
QLOGN = QMIN 

GO. rk: 7 

IF (QLOGN - QMAX) 7, 7, 9 
QLOGN = QMAX 

CONTINUE 

PRINT 120, QLOGN 

FORMAT (10X,4F10.3) 
FORMAT (20X,F10.3) 


END 
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THE TRIANGULAR DISTRIBUTION VARIABLE GENERATOR 


The FORTRAN IV Computer program for deviates with a 


Triangular distribution requires one uniformally distributed 


random function, RANF(O). This uniformally distributed random 


function will assume either the value of W(x < xi), - WXLTXI - 


or W(x > xX; ) - WXGTXI - depending on whether the calculated 


value of x. (XT) 1s equal to or less than the inputed value of b. 


The FORTRAN IV program incorporating these concepts is 


listed as follows: 


50 
ro 


Ha Z 
300 


PROGRAM TRIA (INPUT,OUTPUT) 


THIS PROGRAM CALCULATES THE BEST ESTIMATE OF PROFIT/ BBE 
ASSUMING A * TRIANGULAR * DISTRIBUTION OF PROBABILITY. 


A = MINIMUM PROFIT/BBL EXPECTED 
B = MOST LIKELY VALUE OF PROFIT/BBL EXPECTED 
C = MAXIMUM PROFIT/BBL EXPECTED 


READ: tik 2eec eB nC) 

CALL RANF(0) 

WXLTXI = RANF(O) 

WXGTXI = 1.0 - WXLTXI 

XI = A + SQRT(WXLTXI* (B-A)*(C-A)) 
ie (Ow) ES St 

XI = C - SQRT(WXGTXI* (C-B)*C-A)) 
PROFIT = XI/100. 

PRINT 300, PROFIT 

FORMAT (10X,3F10.2) 

FORMAT (15X,F4.2) 

END 
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THE BETA DISTRIBUTION VARIABLE GENERATOR 


The FORTRAN IV Computer program for deviates with a 


Beta distribution is listed as follows: 


0) 
500 


10 
20 


50 


PROGRAM BETA (INPUT ,OUTPUT) 


THIS PROGRAM CALCULATES THE BEST ESTIMATE OF PROFIT/BBL 
ASSUMING A * BETA * DISTRIBUTION OF PROBABILITY. 


A = MINIMUM PROFIT/BBL EXPECTED 
B = MOST LIKELY VALUE OF PROFIT/BBL EXPECTED 
C = MAXIMUM PROFIT/BBL EXPECTED 


RED Oe) ae). 1c 

XMU = (4.*B+A+C)/6.0 

XVAR = (C-A)**2/36.0 

BMEAN = (XMU-A)/(C-A) 

BVAR = XVAR/(C-A)**2 

XK1 = BMEAN* (BMEAN* (1.0-BMEAN)/BVAR-1.0) 
XK2 = XK1*((1.0-BMEAN) /BMEAN) 

CALL GAMMARN (XK1,GAM1) 

CALL GAMMARN (XK2,GAM2) 

BETA = (GAM1/(GAM1 + GAM2))*(C-A) +A 
PROFIT = BETA/100. 

PRINT 300, PROFIT 

FORMAT (10X,3F10.2) 

FORMAT (15X,F4.2) 

END 


SUBROUTINE GAMMARN (TK, GAM) 
GAMMA = 1.0 

Kl = TK 

TK1l = Kl 

Yl = RANF(O) 

Ter ean Lee GD) ae tn OZ 0 
ie = ia 

BNO! SO Nea al 

Y2 = RANF(0) 

GAMMA = GAMMA*Y2 

GAM = -ALOG (GAMMA ) 

RETURN 

END 





AEPENDIX BD 


COMPUTER PRINTOUT FOR BASE DATA 





PROFIT DISTRIBUTION OF AN EXPLORATORY DRILLING PROGRAM 


NO. OF WELLS PROR OF ORY HOLE 
DRILLED SUCCESS cosT QMEDIAN Q84 QMAX QMIN 
20 e150 e150 20-00 46.00 140.00 1.00 


THE TRIANGULAR DISTRIBUTION FOR DETERMINATION OF THE 
PROFIT PER BARREL HAS THE FOLLOWING EXTREMES IN CENTS = 
A = 50.00 B = 88.00 C = 115.00 


ITERATION NUMBER OF FIELD PROFIT PER’ ITERATION BEST ESTIMATE 


SUCCESSES SHZES BARREL PROFIT OF PROFIT 
l 2 52.90 083 41.17 41.17 
2 3 64.03 1.01 61.90 Sh. 53 
3 0 0-00 0-00 -3.00 33.36 
& & 114.52 le2 114.13 53-55 
5 5 189.57 079 148.40 72-52 
6 0 0-00 0-00 -3-00 59.93 
7 3 34.212 087 27.05 55-24 
8 5 254.88 1-651 254.59 80-15 
9 3 198.18 08] 157.93 88.80 
10 2 32-61 1-01 30-36 82.95 
11 0 0.00 0-00 -3.00 15214 
le 2 38.72 093 33-17 71264 
13 5 179.08 093 164.57 78.79 
14 3 95.99 057 51.88 76.87 
15 3 166.86 61 98.92 78.34 
16 3 129.24 °65 81.67 78.55 
17 2 28.20 86 21.46 13. YS 
18 3 54.65 e80 41.00 73229 
19 3 713-73 095 67.13 72.97 
20 ] 85.62 97 80-58 a, 35 
el 5 127.4] e81 100.96 74.66 
22 1 3284 92 «69 71-30 
23 2 4045 01/9 29-10 69.46 
24 4 59.44 087 49.05 68.61 
25> 2 53-91 081 41.17 67.52 
26 3 49.27 °66 30-22 66.08 
el 2 26-85 e/0 16-03 64.23 
28 5 96.30 ef3 67294 64.36 
29 2 27.60 074 17-80 62.75 
30 + 187.40 099 183.55 66./8 
31 2 137.95 064 85.37 67,38 
32 4 83.15 82 65./0 67.33 
33 2 107.08 82 84.63 67.85 
34 0 0.00 0.00 -3.00 65.77 
35 4 84.91 °85 69.90 65.89 
36 2 38.40 094 33-52 64.99 
37 2 25244 1.92 23.35 63.86 
38 & 125.78 e73 89.62 64.54 
39 l 6.76 1/8 2241 62.95 
40 5 110.73 89 95.83 63.77 
4] 0 0.00 0.00 -3.00 62.14 
42 4 


116.96 1.09 124.90 63.63 
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229.75 
32.97 
25.19 

0.00 

181.65 
84676 

138.59 

160.84 
91.12 
50-01 
30-00 
4341 

129.75 
19.93 
6336 
72.45 

110.63 

100.93 

194.06 

170.74 

0.00 
74.02 
39.45 
89.92 

112.67 

172622 
5186 
4719 
63.76 
76.24 
23-86 

110.39 
20276 
95.34 
38.74 
23.97 
32.79 
S087 
9119 
96.5] 

1i9ea0 

0.00 

18412 
31.18 
51-86 
80.81 

109.49 
S549 
19.07 
29 024 
88.87 
23-84 
83.49 
88.52 

156023 

0-00 
50-68 
42 256 

0-00 
55247 
88.88 
34.29 


09] 


1.07 
0.00 
1-02 

087 


059 
1-04 


1-04 
e7l 
067 
055 
094 


079° 


92 
087 
e63 
e66 
094 
70 
0-00 
e 86 
087 
96 
°84 
88 
e99 
084 
091 
058 
80 
80 
084 
093 
°87 
© 76 
e779 
89 
e73 
°63 
091 
0.00 
e733 
0/8 
86 
058 
e95 
eS/ 
Lele 
094 
e/0 
88 
067 
1-08 
84 
0-00 
1.05 
093 
0-00 
087 
9] 
e89 


206-55 
29-02 
24220 
-~3.00 

183-29 
10292 
79.57 

165.66 
92252 
32264 
17-45 
20-90 

119.35 
12.84 
55-84 
60.17 
66.80 
64.42 

181.04 

117.80 
-3.00 
60-82 
31-43 
83.53 
91.91 

150.10 
48.86 
36-78 
55-78 
41.86 
16.45 
85.40 
14.53 
85.9] 
30-88 
15.85 
23-31 
42.19 
64.26 
58.39 

106.21 
-~3.00 

132.41 
21.46 
41.95 
44.18 

101-59 
28.98 
18.70 
24286 
59.78 
18.05 
53-13 
92.78 

129.02 
-3.00 
50.79 
37.08 
-3.00 
45.67 
18-40 
27059 


66.96 
66.10 
65.17 
63.68 
66.23 
66.33 
66.60 
68.58 
69.05 
68.35 
67.39 
66.53 
67.49 
66.5] 
66.32 
66.22 
66.23 
66.20 
68.08 
68.88 
67.74 
67.63 
67.08 
67.32 
67.69 
68.90 
68.61 
68.16 
67.98 
67.62 
66.92 
67.17 
66.47 
66.72 
66.26 
65.61 
65.08 
64.80 
64.79 
64.71 
65.21 
64.40 
65.20 
64.69 
64.43 
64.20 
64.62 
64.23 
63.72 
63.30 
63.26 
62.78 
62.68 
63.00 
63.68 
63.00 
62.87 
62.61 
61.96 
61.81 
61.97 
61.64 
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105 
106 
107 
108 
109 
110 
111 
112 
113 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
130 
131 
132 
133 
134 
135 
136 
137 
138 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
15] 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
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124.19 
82-08 
92.60 

141.82 

101.35 
45226 
16250 
27el2 
18.92 
83.56 
31.66 
85 46 

166.72 

199.84 

118.33 

150.97 

0-00 

C20. >> 

199.11 

111.75 
25256 

0.00 
57.289 
83.449 

245-33 

194.04 

108.37 
2.57 
18.69 
59.56 

200-74 

128.45 

3-87 

231-8 
46-63 
31.68 
$2013 

100-82 
21243 

0-00 
860 

24661) 
39053 

209-6) 
84.89 
3316 
37.88 

5058 

266254 
47.87 
3016 

196.48 

161.99 
89-22 
40-85 

135.26 
31.66 
95.17 
19.270 

5240 
90-85 
66-02 


e68 
1e13 
e65 
efS 
087 
ofl 
off 
°99 
°96 
83 
e 84 
064 
° 84 
68 
094 
°85 


60 
093 
2 84 
e80 
0.00 
e73 
1.93 
095 
e95 
1-02 
054 
lel 
e86 
092 
095 
«88 
074 
067 
66 
«88 
1-06 
1-10 
0-00 
079 
093 
ef} 
099 
e66 
e775 
e85 
°65 
e7) 
0 f4 
°e%6 


79 


© 83 
e779 
e80 
1-06 
09} 
8% 
99 
e85 
e73 


82.68 
90-48 
57292 
103.75 
85-45 
32-31) 
56434 
24204 
15.36 
67-61 
24.00 
52.203 
137.93 
133.63 
108-42 
125-42 
-3.-00 
151.98 
184.232 
91.60 
17.70 
-3.00 
39242 
83.96 
229278 
181.68 
108.45 
36264 
16.22 
48.53 
182.97 
119.53 
225 
169.64 
28-36 
18.03 
43.28 
104.06 
27-35 
-3.00 
3-96 
227-67 
25-60 
205-16 
53-98 
22037 
29-71 
280 
187-86 
32.89 
26-68 
1S$2.4]) 
159.09 
731230 
29-58 
105.98 
30-89 
84.64 
64.64 
2250 
14.85 
4Se 18 


61.84 
62211 
62.06 
62245 
62.266 
62.38 
62.33 
61.99 
61.58 
61.63 
61.30 
61.22 
61.88 
62249 
62.87 
63-39 
62.84 
63.57 
64.56 
64.77 
64.40 
63.86 
63-67 
63.83 
65211 
66.01 
66.34 
66.11 
65.74 
65.6] 
66.48 
66.87 
66.38 
67.13 
66.85 
66.50 
66.234 
66.60 
66.33 
65.85 
65-42 
66653 
66.25 
67.19 
67.10 
66.81 
66.56 
66.13 
66.92 
66.70 
66.44 
67.00 
67.58 
67.61 
67.37 
67.61 
67.38 
67.49 
67647 
67.07 
67.12 
66.99 
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167 
168 
169 
170 
171 
17/2 
173 
174 
175 
176 
177 
178 
179 
180 
18] 
182 
183 
184 
185 
186 
187 
188 
189 
190 
19] 
192 
193 
194 
195 
196 
197 
198 
199 
200 
201 
202 
203 
204 
205 
206 
207 
208 
209 
210 
211 
212 
213 
214 
2 
216 
2l7 
218 
219 
220 
22) 
222 
223 
224 
225 
226 
2el 
228 
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84.28 
87.96 
48.4) 
47.52 
4962 
9974 
0.00 
57 047 
163-81 
15576 
18530 
0-00 
0-00 
108.78 
82.87 
68.33 
43.93 
92.28 
49077 
10-66 
101-08 
8432 
116.27 
46074 
6346 
69.48 
33-12 
0.00 
127.19 
146.63 
48.18 
81.22 
76037 
4629 
49299 
56042 
1027 
33-89 
12.06 
6732 
146.58 
78.10 
221.96 
32.79 
179.39 
102.93 
182.99 
2319 
179.31 
40.05 
68.35 
94.78 
171.29 
0.00 
110.02 
42210 
142.72 
132.25 
60 006 
1624 
6540 
72.30 


1.05 
89 
e96 
°95 
° 80 
e92 

0-00 
°99 
07S 

1-00 
°68 

0-00 

0-00 
097 
78 
81 
099 
°62 
ee 
82 
¢80 

1-00 
094 
ovf 
78 
e775 
86 


°85 
094 
e790 


078 
e/6 
054 
e/0 
81 
°89 
78 
087 
©90 
° 70 
e7l 
78 
°62 
°65 
059 
°90 
°85 
°89 
e179 
90 
097 


°92 
68 
e95 
e82 
70 
e/0 
094 
1.04 


8629 
75-31 
44.218 
42.63 
37204 
89.48 
~3.00 
54-24 
120-32 
153.48 
124.21 
~3.00 
~3.00 
102.99 
61.96 
52.239 
40.75 
54.36 
32.94 
592 
78-87 
81.58 
10632 
33-63 
47-11 
49-50 
25/2 
~3.00 
105.54 
135.76 
3124 
81.18 
57-05 
32242 
24.2) 
36-53 
5-52 
27-39 
6259 
55-9e 
128-85 
51-99 
155-14 
23-02 
109.95 
65.00 
106-45 
18-16 
149.58 
33-09 
51-88 
82.79 
164.52 
~3.00 
98.94 
26-12 
132-74 
105.64 
39254 
8.69 
58.74 
72230 


67.210 
67.15 
67.202 
66.87 
66.70 
66.83 
66.43 
66.36 
66.66 
67.16 
67.48 
67.08 
66.69 
66.89 
66.87 
66.79 
66.65 
66.58 
66.40 
66.07 
66.14 
66622 
66.43 
66.26 
66.216 
66-07 
65.87 
65.51 
65.72 
66.07 
65.90 
65.97 
65.93 
65.76 
65.55 
65-41 
65ele 
64.93 
64.65 
64.60 
64.9] 
64.85 
65-28 
65.08 
65.30 
6529 
65-49 
65.27 
65.66 
65.51 
65.44 
65.52 
65.98 
65.66 
65.81 
65.63 
65.294 
66.11 
65.99 
65.74 
65.71 
65.274 
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230 
231 
232 
233 
234 
235 
236 
237 
238 
239 
240 
241 
242 
243 
244 
245 
246 
247 
248 
249 
250 
251 
252 
253 
254 
255 
256 
257 
258 
259 
260 
261 
262 
263 
264 
265 
266 
267 
268 
269 
270 
271 
272 
273 
274 
275 
276 
277 
278 
279 
280 
281 
282 
283 
284 
285 
286 
287 
288 
289 
290 
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72.06 
128.08 
4547 
200.29 
38.19 
67.98 
43240 
121.19 
30296 
70.22 
90 32 
0.00 
84.62 
63.96 
119.11 
75.18 
163.88 
62.32 
43.48 
82.25 
51.07 
35.34 
62.71 
15.34 
113.74 
62.76 
183.22 
51.97 
32 48 
95.57 
24.58 
100.22 
126.84 
21242 
14.90 
81.39 
192.19 
0.00 
166.32 
115.73 
237611 
82.53 
27.87 
28230 
10.30 
136.50 
97233 
110.15 
117.57 
56.03 
66687 
104.59 
81.22 
56046 
102.32 
0.00 
135.46 
150.39 
20-62 
25.73 
39.28 


78 
°89 
67 
1-01 
75 
«86 
7/6 
84 
e73 
61 
e66 
0.00 
°92 
82 
063 
e910 
e85 
e 86 
80 
96 
of4% 
059 
98 
96 
80 
°62 
76 
1-658 
ef/0 
1-01 
e 86 
87 
096 
95 
285 
90 
1-04 
0.00 
075 
e 76 
098 
093 
286 
ell 
o77 
084 
«86 
ell 
073 
e/0 
ofS 
e65 
of4 
1-10 
1-10 
0.00 
«80 
0ST 
°65 
1.05 
69 
1.05 


53-48 
111.51 
27-83 
200.10 
25-98 
56-05 
30-58 
99.09 
19.73 
40.32 
57-49 
-3.00 
15216 
49.89 
T2276 
64.79 
136.17 
51.07 
32-27 
76280 
35.23 
18.11 
59-10 
11.95 
88.76 


36.200 
137.64 


53-87 
19.95 
93.89 
18.21 
85.22 
119.65 
17.54 
9.85 
71.00 
197.50 
-3.00 
121.245 
85.21 
230-96 
73-67 
21-19 
19.11 
5-10 
111.76 
80-92 
83.08 
83.53 
36.279 
47.18 
65.00 
57.78 
59.64 
109.98 
-3.00 
105.93 
84.00 
10.98 
176.30 
15.04 
38.254 


65.69 
65.89 
65.72 
66.30 
66.13 
66.08 
65.93 
66.07 
65.88 
65.77 
65.74 
65.45 
65.49 
65.42 
65.45 
65.45 
65.74 
65.68 
65.55 
65.59 
65.47 
65.28 
65.26 
65.04 
65.14 
65.02 
65.31 
65.26 
65.09 
65.20 
65.02 
65.09 
65.39 
65-12 
64.91 
64.93 
65.43 
65.18 
65.39 
65.46 
66.08 
66.11 
65.94 
65.77 
65.55 
65.71 
65.77 
65.83 
65.90 
65.79 
65./2 
65./2 
65.69 
65.67 
65.83 
65.59 
65.73 
65.79 


65.60 


65.98 
65.81 
65.71 
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291 
292 
293 
294 
295 
296 
297 
298 
299 
300 
301 
302 
303 
304 
305 
306 
307 
308 
309 
310 
311 
312 
313 
314 
315 
316 
317 
318 
319 
320 
321 
322 
323 
324 
325 
326 
327 
328 
329 
330 
331 
332 
333 
334 
335 
336 
337 
338 
339 
340 
341 
342 
343 
344 
345 
346 
347 
348 
349 
350 
351 
352 
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111.95 
46.71 
148.66 
42.49 
72 296 
58.64 
65.29 
59.52 
4119 
45.93 
111235 
81-95 
9202 
25 266 
208 044 
178.62 
58.22 
305.88 
318.55 
S894 
9.99 
144.74 
290035 
90-82 
90-02 
46678 
66-56 
712) 
35.39 
56.28 
21659 
77244 
12-08 
111-9) 
33-23 
130-08 
0-00 
36-83 
198.73 
127-67 
8149 
10-50 
135.69 
110.81 
S0.27 
150.83 
109.13 
10749 
52.73 
116.25 
2424) 
4724 
164.77 
165.69 
210.82 
140.86 
118.45 
39.52 
13567 
45.99 
56-88 
53.83 


e775 
090 
90 
054 
°63 
°75 
e91 
e95 
°92 
1-01 
064 
°82 
09] 
093 
059 
e85 
89 
e84 
88 


087 
81 
88 
89 
°95 
059 
e7f 
093 
efS 
e89 
094 
093 
ofl 
«82 
e 88 
093 
0-00 
0/8 
e92 
°89 
°95 
98 
81 
1-02 
80 
e88 
88 
e82 
e62 
lele 
81 
«80 
©96 
e 84 
e 80 
°66 
057 


89 
°85 
073 
79 


81.07 
39259 
130.68 
20228 
43-69 
41/2 
56 88 
54el2 
35-28 
43-62 
68244 
64.63 
5-48 
2lel? 
121.06 
149.63 
49.17 
254.27 
279-47 
56-88 
6.03 
114.40 
252-48 
78280 
82-82 
25014 
49.€8 
64.633 
23-6) 
47.70 
1749 
69272 
6-66 
88.91 
26246 
117-92 
-3.00 
26-00 
179.64 
112.09 
75.209 
7239 
107-22 
109.48 
37-68 
129.77 
93.81 
86.15 
30-55 
128.4) 
16.91 
35245 
156.25 
137.206 
166.73 
90.53 
65.08 
38-94 
117.93 
36244 
39231 
39272 


65.77 
65.68 
65.90 
65.74 
65.67 
65.59 
65.56 
65.52 
65.42 
65.35 
65.36 
65.35 
65.16 
65.01 
65.20 
65.47 
65.42 
66.03 
66.72 
66.69 
66.50 
66.65 
67224 
67.28 
67.33 
6720 
6le 14 
67.13 
66.99 
66.93 
66.78 
66.79 
66.60 
66.67 
66-55 
66-70 
66.49 
66.37 
66.71 
66.85 
66.87 
66.69 
66.82 
66.94 
66.86 
67.04 
67.12 
67.18 
67.07 
67225 
67.10 
67.01 
67.27 
67.47 
67.76 
67.83 
67.82 
67.74 
67.88 
67.79 
67371 
67.63 
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354 
355 
356 
357 
358 
359 
360 
361 
362 
363 
364 
365 
366 
367 
368 
369 
370 
371 
372 
373 
374 
375 
376 
377 
378 
379 
380 
381 
382 
383 
384 
385 
386 
387 
388 
389 
390 
391 
392 
393 
394 
395 
396 
397 
398 
399 
400 
401 
402 
403 
404 
405 
406 
407 
408 
409 
410 
411 
412 
413 
414 
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36215 
85.85 
63240 
23-07 
24240 
233-50 
40-20 
99 204% 
10-86 
42047 
140.00 
31-09 
160.49 
11.79 
116.60 
114.60 
19.236 
148-53 
118.97 
15.239 
123.28 
229-61 
89233 
66.18 
19-08 
40-68 
Se2l 
72-81} 
50-20 
30-01 
6.09 
52-00 
67-27 
1025 
97-62 
110.19 
96.36 
10-42 
125.93 
60-80 
26.17 
71.205 
139.19 
10-88 
157.85 
7256 
56-47 
187.61 
93-10 
7209 
21.243 
106.234 
67.93 
61.13 
106.43 
118.99 
14.24 
267.210 
68.07 
34.285 
0-00 
7-86 


097 
73 
81 
°/6 
1-04 
1.02 
e 86 
ef9 
ele 
91 
82 
89 
°85 
°99 
e85 
84 


1-08 
1.00 
1.00 
°92 
1.95 
°85 
°/8 
© 84 
e82 
°82 
off 
°66 
1-04 
1.03 
°99 
084 
°95 
080 
84 
087 
064 
ef3 
80 
ef2 
8) 
89 
e5] 
92 
e80 
087 
of4 
88 
90 
68 
1-01 
1.03 
067 
88 
°85 
1.10 
053 
89 
°90 
0.00 
82 


3244 
60.06 
48.92 
14-67 
22.73 
235-80 
31275 
1S e22 
4.93 
35293 
111.65 
25-06 
133.43 
8.79 
96.71 
93243 
13-08 
158.49 
116.98 
72285 
111.40 
239213 
V3e77 
48-92 
13.236 
30-75 
1245 
53-90 
30-298 
28424 
3240 
48.82 
S46 11 
6-87 
14.98 
90.93 
81-36 
42234 
89-71 
46.07 
15.86 
55.27 
121.76 
22/0 
142.81 
3023 
46.54 
136.74 
19202 
3254 
11-81 
104.92 
67244 
38.03 
91.25 
98.09 
12.89 
138.232 
58.28 
28-44 
-3.00 
3-58 


67.53 
67.51 
67.246 
67.31 
67.18 
67.65 
67.55 
67.58 
67.240 
67.32 
67644 
67.232 
67.50 
67.34 
67242 
67649 
67.235 
67.59 
67272 
67.74 
67.86 
68.31 
68.33 
68.28 
68.13 
68.03 
67.86 
67.82 
67.72 
67.62 
67.45 
67.40 
67.37 
67.2) 
67623 
67.29 
67.33 
67.26 
67.232 
67.2ec7 
67214 
67.11 
67.24 
67.08 
67.27 
67.11 
67.206 
67.223 
67.26 
67.11 
66.97 
67.06 
67.06 
66.99 
67.05 
67.213 
66.99 
67-17 
67.15 
67.05 
66.88 
66.73 
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415 
416 
417 
418 
419 
420 
421 
422 
423 
424 
425 
426 
427 
428 
429 
430 
43) 
432 
433 
434 
435 
436 
437 
438 
439 
440 
44} 
442 
443 
444 
445 
446 
447 
448 
449 
450 
45) 
452 
453 
454 
455 
456 
457 
458 
459 
460 
461 
462 
463 
464 
465 
466 
467 
468 
469 
470 
471 
472 
473 
474 
475 
476 
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22.39 
103-76 
69.77 
132.76 
31.52 
15271 
44.28 
54.03 
74.81 
11.22 
73613 
50.83 
22.22 
79.78 
42.77 
189.88 
98.89 
199.46 
187.53 

0-00 

6-00 
62.52 
89.47 

9-12 
120.28 
7350 
115015 
35.23 
34.12 
62.27 
63.79 
95.88 
77292 
83.44 
178.90 
37.46 
113-76 
35021 
34.86 

0-00 
3649 
29.00 
42.32 
45.82 
33.20 
20.75 
Wit 
94.48 
19.87 
44.26 
52.04 
66.38 
110.36 
SS .83 
28.97 
53.71 
62.70 
89.33 
109.69 
45.57 
184.99 
173.40 


1-03 
e86 
6) 
0S7 
097 
090 
094 
90 
95 

Lefe 
°92 
e932 
095 
e80 
e89 
90 
058 
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796 
wi 
798 
799 
800 
80] 
802 
803 
804 
805 
806 
807 
808 
809 
810 
811 
812 
813 
814 
815 
816 
817 
818 
819 
820 
821 
822 


823 


825 
826 
827 
828 
829 
830 
831 
832 
833 
834 
835 
836 
837 
838 
839 
840 
841 
B42 
843 
B44 
845 
B46 
847 
848 
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268 235 
52216 
2212-35 
123.37 
220695 
17259 
0.00 
62.05 
149.27 
64255 
196.83 
3987 
121.40 
81.234 
0.00 
64.46 
35-98 
3220 
59.10 
107.240 
100-10 
42-48 
22294 
143-80 
39639 
121.03 
7.69 
l02.77 
150.14 
53-61 
101.02 
129.17 
41.61 
103.07 
11240 
135-28 
148.40 
45-92 
147220 
160-11 
310218 
B2ei7 
115.43 
10.69 
143-61 
121-62 
141.98 
153.20 
0-00 
214.8) 
16.05 
80 43 
26.99 
35.292 
119.20 
17649 
99.31 
15.99 
133.08 
80.97 
16.01 
75.06 


069 
e86 
058 
e/70 
°99 
88 
0-00 
085 
e67 
ef2 
079 
1-02 
e65 
ae | 
0.00 
0/8 
1.05 
ef2 
2/8 
058 
°83 
0f4 
ele 
098 
o/8 
2 84 
ao? 
e 86 
e 89 
053 
e91 
1.%4 
e 84 
084 
067 
e86 
084 
1-01 
e 84 
96 
093 
098 
of9 
095 
1-99 
ef4 
82 
°/6 
0.00 
91] 
°86 
e779 
e59 
81 
090 
056 


097 
°85 
e87 
°99 
°65 


183.17 
42.73 
129.15 
83-23 
216-57 
65.76 
-3.00 
49.93 
97-15 
44.205 
154.07 
58-02 
166592 
11.68 
-3.00 
48.03 
3522 
20047 
43.69 
59.73 
80.88 
28.83 
13-76 
139.1% 
28.00 
99.36 
4.03 
86.06 
13044 
25-83 
89.00 
131./¢5 
32013 
84.30 
73-00 
114-57 
123-12 
43-66 
121-28 
152.26 
286-23 
77-82 
89.25 
1229 
154.52 
87214 
113.49 
113.75 
-3.00 
193-06 
11.00 
60-80 
12.99 
26-60 
105.59 
72-15 
103.34 
12.70 
111.29 
68.18 
13-22 


46-46 


67290 
67.86 
67.94 
67.96 
68.15 
68.15 
68.06 
68.03 
68.07 
68.04 
68.15 
68.11 
68.12 
68012 
68.04 
68.01 
67.97 
67.91 
67.88 
67.87 
67.89 
67.84 
67.77 
67.86 
67.81 
67.85 
67.77 
67.79 
67.87 
67.82 
67.84 
67.92 
67.88 
67.90 
67.91 
67.96 
68.03 
68.00 
68.06 
68.17 
68.43 
68.44 
68.47 
68.39 
68.50 
68.52 
68.57 
68.63 
68.54 
68.69 
68.62 
68.61 
68.55 
68.50 
68.54 
68.47 
68.51 
68.44 
68.49 
68.49 
68.43 
68.40 





849 
850 
851 
852 
853 
854 
855 
856 
857 
858 
859 
860 
861 
862 
863 
864 
865 
866 


868 
869 
870 
87) 
B72 
873 
874 
875 
876 
877 
878 
879 
880 
881 
882 
883 
884 
885 
886 
887 
888 
889 
890 
891 
89e 
893 
894 
895 
896 
897 
898 
899 
900 
901 
902 
903 
904 
905 
906 
907 
908 
909 
910 
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16649 
131.13 
55-97 
175-05 
148.64 
155-51 
99.10 
61.57 
45298 
42-86 
57-38 
278.61 
7282 
21-67 
89.49 
132.26 
16-65 
15261 
30 31 
59.10 
51.81 
78.80 
35014 
42.67 
92ele 
25-99 
7283 
98.31 
60-85 
82245 
4622 
105.05 
72259 
131.70 
42-01 
55.18 
49.31 
33-63 
123-45 
$2231 
133-15 
118.14 
117-93 
175.85 
89.84 
14.98 
60-13 
76265 
36-35 
18.17 
60.94 
69.61 
165.45 
44.09 
121-40 
32284 
88.36 
6.56 
299.06 
67.80 
98.84 
64.58 


88 
1.08 
1.02 
085 
064 
275 
1.02 
°68 
°79 
o77 
°73 
80 
71 
ste 
297 
290 
88 
265 
°63 
290 
80 
°82 
85 
289 
269 
1.09 
255 
°79 
kf 
91 
064 
1204 
261 
98 
1-03 
295 
°62 
1607 
84 
1012 
291 
290 
299 
88 
.78 
1204 
°79 
60 
92 
lel] 
294 
96 
6] 
279 
81 
99 
91 
265 
292 
~90 
7 
296 


65.20 
139.69 
54.267 
145.97 
93243 
114.18 
98.70 
39.06 
33-65 
30235 
39 48 
219.98 
2e7) 
12.95 
84.4) 
115.94 
11.73 
464) 
16-38 
50.14 
38.95 
62.07 
27239 
35-67 
61.214 
25.258 
1.49 
T4279 
41.18 
72205 
27-37 
107274 
41.54 
126-04 
40-67 
49.75 
27.65 
33-22 
101.33 
44.13 
119.60 
104-31 
113.85 
153.03 
67-39 
12.70 
45.10 
43-22 
30.89 
17.27 
594.87 
63.87 
98.87 
32-31 
95.86 
29.82 
77291 
1.42 
274.07 
58.48 
67.179 
59.49 


68.40 
68.48 
68.47 
68.56 
68.59 
68.64 
68.67 
68.64 
68.60 
68.55 
68.52 
68.70 
68.62 
68.55 
68.57 
68.63 
68.56 
68.54 
68.48 
68.46 
68.42 
68.41 
68.37 
68.33 
68.32 
68.c7 
68.20 
68.20 
68.17 
68.18 
68.13 
68.18 
68.15 
68.21 
68.18 
68.16 
68.11 
68.07 
68.11 
68.09 
68.14 
68.18 
68.23 
68.33 
68.33 
68.27 
68.24 
68.21 
68.17 
68.11 
68.10 
68.10 
68.13 
68.09 
68.12 
68.08 
68.09 
68.02 
68.24 
68.23 
68.23 
68.22 





911 
912 
913 
914 
915 
916 
917 
918 
919 
920 
921 
922 
923 
924 
925 
926 
927 
928 
929 
930 

931 

932 
933 
934 
935 
936 
937 
938 
939 
940 
941 

942 
943 
944 
945 
946 
947 
948 
949 
950 
951 
952 
953 
954 
955 
956 
957 
958 
959 
960 
961 

962 
963 
964 
965 
966 
967 
968 
969 
970 
971 

972 
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64-88 
15228 
120.56 
103.66 
59el2 
106.89 
117.69 
91.222 
120.53 
71218 
63.55 
41.41 
61.95 
l7e27 
28294 
57-08 
160-68 
44.95 
57239 
285-64 
18-85 
131-95 
24-78 
149.44 
25-00 
108-64 
127-18 
34.01 
68.04 
114.78 
93-84 
12216 
235-57 
141.75 
86-68 
95-48 
12.50 
94.18 
131.53 
46295 
0.00 
41.66 
107.98 
112.16 
58-17 
18-07 
226-90 
71-78 
18.57 
120-67 
T1421 
Gell 
102-85 
48.64 
134.03 
141-04 
58.36 
43.71 
245-90 
118.89 
81.39 
99-31 


1.01 
lell 
095 
79 
90 
095 
98 
096 
68 
69 
097 
e79 
074 
© 64 
1-04 
1-03 
093 
86 
067 
e776 
e 76 
1-01 
1.05 
e82 
° 83 
°8] 
2°60 
1-01 
e176 
83 
074% 
e90 
91 
90 
91 
83 
1-00 
82 
81 
094 
0.00 
1-00 
°69 
°81 
1-00 
73 
°62 
°93 
°81 
1-09 
80 
094 
81 
067 
e88 
1-07 
094 
90 
16 
92 
95 
e6l 


62252 
81.36 
ll2e32 
19.48 
50-52 
98.41 
112.95 
84.89 
78.87 
46.95 
58-82 
3024] 
43.02 
8223 
27233 
56-27 
146.55 
35-79 
36220 
213-59 
11.59 
131-12 
23.11 
119.74 
18.00 
85245 
732-65 
31.54 
49.10 
92-66 
66.8] 
8.12 
212-77 
126.02 
76262 
16282 
9-69 
14629 
104-49 
41.47 
3.00 
39.17 
72-56 
88.83 
55.31 
10.31 
139.19 
64.04 
12-40 
128-42 
56-32 
6246 
80.43 
29-81 
115.00 
148.36 
52-05 
36-96 
184.27 
106-62 
14.41 
63.99 


68.22 
68.23 
68.28 
68.29 
68.c7 
68.30 
68.35 
68.37 
68.38 
68.36 
68.35 
68.31 
68.28 
68.21 
68.17 
68.16 
68.24 
68.2) 
68.17 
68.33 
68.27 
68.34 
68.29 
68.34 
68.29 
68.31 
68.31 
68.27 
68.25 
68.28 
68.28 
68.cl 
68.37 
68.43 
68.44 
68.245 
68.38 
68.39 
68.43 
68.40 
68.32 
68.29 
68.30 
68232 
68.31 
68.25 
68.232 
68.31 
68.26 
68-32 
68.31 
68.24 
68.26 
68.22 
68.26 
68.35 
68.33 
68.30 
68.42 
68-46 
68.46 
68.46 
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973 
974 
975 
976 
977 
978 
979 


981 
982 
983 
984 
985 
986 
987 
988 
989 
990 
991 
992 
993 
994 
995 
996 
997 
998 
999 


1000 
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75264 
6862 
49.27 
120.25 
113-96 
9.93 
23.98 
0.00 
71.260 
36.94 
48.01 
59.59 
91615 
7220 
89.59 
33-57 
50.60 
96.38 
47.68 
91.39 
30218 
170.25 
54044 
237043 
13.36 
61.66 
112.05 
43.11 


1.01 


074 


1-10 
80 
063 

1.07 
074% 
°82 

1.01 
097 
°82 
e77 
me 
99 
058 
90 
e63 


74.07 
64.64 
37295 
114.25 
18262 
62/6 
17.12 
-3.00 
54254 
3093 
29-35 
41.24 
50.75 
5.07 
69.07 
18.46 
51-47 
68.90 
36-88 
89.55 
26677 
138.48 
39.236 
172.56 
10.38 
33207 
98.59 
24454 


68.46 
68.46 
68243 
68.48 
68.49 
68-42 
68.37 
68.30 
68.28 
68.25 
68.21 
68.18 
68.16 
68.10 
68.10 
68.05 
68.03 
68.03 
68.00 
68.02 
67.98 
68.05 
68.02 
68.13 
68.07 
68.03 
68.06 
68.02 


THE BEST ESTIMATE OF PROFIT BY THE MONTE CARLO SIMULATION METHOD 


AT 1000 ITERATIONS 


THE VARIANCE IS 
THE STANDARD DEVIATION 
THE COEFFICIENT OF VARIATION IS 


68.02 MILLION DOLLARS 


28792367 
IS 


53.660 
«7889 





NO. 


PROFIT DISTRIBUTION OF AN EXPLORATORY DRILLING PROGRAM 


OF WELLS 


DRILLED 


20 


THE TRIANGULAR DISTRIBUTION FOR DETERMINATION OF THE 


PROB OF ORY 
COST 
e150 


SUCCESS 
©150 


HOLE 


QMEDIAN 
20-00 


Q84 
46-00 


QMAX 
140.00 


QMIN 
1.00 


PROFIT PER BARREL HAS THE FOLLOWING EXTREMES IN CENTS - 
= 88.00 


A 


50-00 


UPPER LIMIT OF 
PROFIT INTERVAL 


~10 


10 
20 


40 
50 


70 
80 


100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
350 
400 
450 
500 
550 


B 


PROBABILITY 
OF INTERVAL 


#0000 
©9420 
«9530 
©9740 
20910 
- 1060 
00870 
©0780 
©0630 
©0630 
-0650 
©0450 
00490 
09340 
e0280 
0280 
©0130 
e0220 
e0090 
©0050 
0110 
0070 
-0050 
©0050 
©9030 
-0030 
*0000 
©0040 
~0020 
~0030 
©9010 
#0000 
#0000 
20010 
#0000 
#0000 
#0000 


C= 


CUMMULATIVE 


115.00 


PROBABILITY 


#0000 
©0420 
©0950 
~ 1690 
«2600 
« 3660 
©4530 
05310 
05940 
~6570 
e/220 
7670 
«8160 
- 8500 
- 8780 
~9060 
©9190 
©9410 
«9500 
9550 
e 9660 
09730 
©9780 
©9830 
9860 
«9890 
9890 
©9930 
©9950 
09980 
9990 
9990 
9990 
1.0000 
1.0000 
1.90000 
1.0000 
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Lawrence William Vogel was born in Brooklyn, New York 
on January 5, 1938, the son of Charles and Lillian Vogel. He 
attended the Brooklyn Technical High School and the City 
College of the University of the City of New York where he was 
awarded a Bachelor of Civil Engineering degree in June 1959. 

After graduation and until entering OCS in November 
1959, he was employed as a Civil Engineer with the M.W. Kellogg 
Company, New York. He was commissioned Ensign, in April 1960 
and currently holds the rank of Lieutenant Commander, Civil 
Mmesameer Corps, U.S. Navy. 

Subsequent to his commissioning, LCDR Vogel was 
assigned as Project Management Officer and then as Engineering 
Officer at the Public Works Department, Yokosuka, Japan, with 
additional duty in Contract Administration. Following this 
three year assignment, LCDR Vogel was ordered to the Public 
Works Center, Newport, Rhode Island, where he served in a 
facilities management position. For his next assignment, he 
served as Assistant Resident Officer in Chagre of Construction, 
Newport. Following this tour of duty he served as Staff Civil 
Engineer, U.S. Naval Air Station Norfolk, Virginia and in June 
1967 was ordered to the Construction Directorate, U.S. Military 
Assistance Command, Vietnam. While serving with MACV, LCDR 
Vogel was assigned to the Base Development Division where he 
was responsible for construction planning, determining const- 


ruction requirements and for monitoring all construction in 
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the I Corps Tactical Zone tor allerrvec World wiiiiicas) Ss owecer 
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